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SY NOPSIS 

CLASSICAL AND SEMI CLASS I CAL MICROSCOPIC APPROACHES 
FOR HEAVY-ION FUSION REACTIONS 

SUBODH S. GODRE 
Ph. D. 

Department Of Physics 
Indian Institute of Technology, Kanpur 

August, 1988 

In heavy-ion collisions nuclei are subjected to much stronger 
srturbations than it is possible with the lighter projectiles, 
lich opens up channels for new reaction processes. At 
rojectile energies of the order of a few MeV/nucleon, reaction 
rocesses such as deep-inelastic scattering and fusion dominate, 
hich involve large scale rearrangement of nucleons and transfer 
f energy from the relative motion to the internal excitations, 
n this thesis a study of heavy-ion fusion reactions in the 
lassical-microscopic approach and the formulation of a 
emiclassical-microscopic approach are presented. The thesis 
onsists of six chapters. Problems and motivations for the 
resent work are discussed in the first chapter along with an 
verview of the heavy-ion reactions and a brief account of the 
arious theoretical approaches. 

In the classical-microscopic approach nucleons are 
onsidered as classical point particles interacting via suitable 
wo-body forces. The ground-states of the individual nuclei used 
n the collision calculations are first obtained by either 
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proper choice of "the parameters of the NN potential, calculated 
fusion cross sections for 16 0+ 16 0 and 40 Ca+ 40 Ca reactions agree 
well with the experiments. At incident energies above a certain 
threshold, Time-Dependent Hartree-Fock ( TDHF ) calculations 
predict the existence of an 1-window for heavy— ion fusion 
reactions - for which there is no conclusive experimental 
evidence so far. It is shown in chapter 4 that 1-window for 
fusion does not exist in the classical-microscopic approach. 
Calculations in this approach at still higher energies show the 
experimentally observed "side-splash" of nuclear matter in heavy- 
ion collisions - an effect that is completely absent in the TDHF 
calculations due to the neglect of two-body collisions. 

Thus, we find that the classical-microscopic approach 
provides a good description of the heavy-ion fusion reactions, 
without requiring any additional assumptions such as those of 
frictional forces or shape parametrizations as in the case ' of 
macroscopic approaches. However, due to the explicit neglect of 
the zero-point energy we find that the effective NN potential has 
large repulsive-core, smaller depth and narrow width, which gives 
higher values of the compressibility coefficients. Also, due to 
the description of nucleons as point particles, the nuclei do not 
have continuous density distributions. 

On the other hand, we know from the experiments that 
lucleons do have finite size which is not negligible compared to 
the average internucleon separation in the ground-states. 
Therefore, in chapter 5, we replace point nucleons by gaussian 



cyclically minimizing the total potential energy of an initially - 
random distribution of nucleons with respect to small spatial 
displacements of individual nucleon coordinates or by a dynamic 
cooling method. The zero-point motion of nucleons in the ground- 
state configuration is explicitly neglected due to the problem 
of premature evaporation of nucleons from the nuclei. However, an 
effective nucleon-nucleon (NN) potential is used which correctly 
reproduces the ground-state binding energy and the rms radius. 
With the given initial conditions for collisions, trajectories of 
all the constituent nucleons are determined by numerically 
integrating coupled Newton’s equations of motion. Nucleon 
trajectories thus obtained are then analyzed to find fusion cross 
sections. Details of this approach are given in chapter 2. 

Fusion cross sections for some heavy-ion reactions, earlier 
calculated in this approach, are overestimated compared to the 
experimental data at lower energies. Taking ■ LD 0+ °0 reaction as 
an example, these calculations are re-examined in chapter 3. We. 
find that the above discrepancy was due to certain assumptions 
earlier made in the specification of the initial conditions; in 
particular due to the neglect of the Coulomb interaction between 
the colliding ions at far off distances. The use of Lennard- 
Jones form of the NN potential in the earlier calculations is 
also critically examined. 

In chapter 4, a systematic study of the dependence of fusion 
cross sections on the range of a soft-core gaussian form of the 
NN potential shows that fusion cross sections are sensitive to 
j-i-- i — »+•.+. Tflnt. i vfi-tail of the NN potential. With a 
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wave-packets with "frozen" widths and find that the incorporation 
of the finite size of nucleons gives a smoother nuclear density- 
distribution and smaller values of the compressibility 
coefficients compared to the point particle calculations. rp ke 
effective NN potential has smaller repulsive-core and longer 
attractive-tail. However, due to the longer tails, calculated 
values of fusion cross sections for the 16 0+ 16 0 reaction are 
enhanced. 

Finally, in chapter 6, we describe a semiclassical- 

microscopic approach in which it is assumed that the wave-p ackets 
always remain as gaussians, but the "frozen-width" assumption is 
removed. The width of every individual wave-packet is determined 
from the instantaneous configuration of the many wave-p acket 
system by equating the zero-point energy of the wave-packets to 
the Thomas-Fermi kinetic energy in the local-density 

approximation. As a simple example we study the dynamics 0 f a 
two wave-packet system which demonstrates the effects of the 
"Pauli-repulsion" . We also find the ground-state of a m any- k ody 
system. While reproducing the ground-state binding energy an( j 
the rms radius of by adjusting the potential parameters, we 
find that the average kinetic energy per nucleon in the ground- 
state is also in close agreement with the experimental Value 
Thus, the effects of the "zero-point" energy and the . "Pauli- 
repulsion" in the ground-state of the nuclei as well as in 
dynamics of the wave-packets are explicitly incorporated in this 
approach, which can be used for the heavy-ion c °Hision 


calculations . 



CHAPTER 1 


INTRODUCTION 


1.1 HEAVY -I ON REACTIONS AS A NUCLEAR PROBE 

The atomic nucleus was first discovered in 1911 in 

experiments conducted by Rutherford and his associates on 

scattering of a-particles by atoms. Since then for a long time 

till seventies, much of the information on the structure and 

properties of nuclei was obtained by using light projectiles such 
3 

as e,p,n,d,t, He,ot-particles, and pions as probes with the 

reaction mechanisms such as elastic and inelastic scattering, 

compound nucleus reactions and direct reactions such as 

pick-up, break-up, stripping, charge -exchange, few nucleon transfer 

etc.* Only in the early seventies a variety of projectiles with 

heavier masses and a wide range of energies became available with 

2 

the advent of heavy-ion accelerators. 

The main motivation for carrying out experiments with heavy- 

3 

ion beams, however, was a desire to produce superheavy-elements. 
The existence of superheavy-elements was predicted in 1966 by 
Myers and Swiatecki 4 using a liquid drop mass formula and the 
empirical shell corrections to the nuclear binding energy. Based 
on this model it was predicted that significant fission barriers 
exist for nuclei that lie near the magic numbers, Z=114 and N=184. 
It was felt that the best way to produce superheavy- elements would 
be to try to fuse two heavy nuclei under such conditions that the 



2 


fused system would not break apart into pieces. However, all the 
experimental efforts to search for superheavy-elements in heavy- 
ion reactions have failed so far . 5 

Though, the discovery of superheavy-elements has still 

remained elusive, the availability of heavy-ion projectiles has 

opened up new frontiers in the study of nuclear matter. In 

heavy-ion collisions, nuclei are now subjected to much more severe 

perturbations than were earlier possible with the lighter 

projectiles. Heavy-ion projectiles colliding with light and 

heavy target nuclei bring along with them high values of charge, 

mass, kinetic energy and angular momentum, thus opening up 

channels for new reaction processes involving high charge and mass 

transfer, deep-inelastic collisions, or formation of a highly 

6 _ 8 

excited and rapidly spinning compound nucleus. 

So far it has not been possible to reach the predicted island 

of stability around the superheavy-elements in the periodic table. 

Nevertheless, with the help of various reactions involving heavy 

ions, some trans-uranium nuclei and a large number of 

neutron-rich and proton-rich nuclei, far away from the line of 

beta- stability, have been produced, thus probing the limits of 

nuclear stability . 10 Involvement of large angular momenta in 

heavy-ion collisions is helpful in studying the limit of stability 

11 12 

of nuclei under fast rotations and strong deformations. ’ In 

heavy-ion collisions at intermediate and high projectile energies, 
high densities and temperatures of nuclear matter are 
reached . 13 ’ 14 The behaviour of nuclear matter at varying 
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densities and temperatures can throw some light on the 
equation-of -state of nuclear matter 13,14 which, besides its 
intrinsic interest is also essential for the understanding of the 
properties of the astrophysical objects such as neutron stars. At 
very high relativistic energies one can form a hot superdense 
piece of nuclear matter with the possibility of a phase transition 
of the baryonic matter to a new phase of matter called the 
Quark-Gluon-Plasma. 

Thus , with the help of the new tools in the form of heavy 

ions at disposal, one can now reach those regions of nuclear 

physics which were not accessible in the past with the lighter 

projectiles. Apart from the domain of Nuclear Physics, heavy-ion 

collisions are also being used for studies in the domain of 

Atomic Physics, and also have important practical applications in 

ion-implantation techniques, semiconductor doping, radiation 

2 

biology and therapy etc. 

Since heavy-ion systems have only few tens to few hundreds of 
nucleons depending on the target-projectile combinations, these 
systems are neither simple few-body systems nor many-body systems 
in the sense of a gas. Thus, theoretical understanding of the 
heavy-ion reaction phenomena also presents a challenging task. In 
section 1.2 we give a classification of heavy-ion reactions. Some 
of the approaches in studying heavy-ion collisions and motivations 
for the present study are discussed in section 1.3. 
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1-2 CLASSIFICATION OF HEAVY -I ON COLLISIONS 

A general physical description and the classification of the 

many processes occuring in heavy-ion collisions is given below. 

Overall features of heavy-ion collisions can be described by 

1 R 

reference to Figure 1.1. 

At collision energies below the Coulomb-barrier the ions do 
not come in contact and interact only through the Coulomb field, 
resulting in Rutherford scattering and possibly Coulomb excitation 
of the target and projectile nuclei. Quantum mechanically some 

1 Q 

tunneling may take place leading to sub-barrier fusion. 

At energies above the Coulomb-barrier the ions can come 
closer and feel the strong nuclear potential. It then becomes 
more convenient to discuss the interaction in terms of the impact 
parameter (b) and one can broadly classify different reactions as 

follows : 

(1) If the impact parameter b is larger than the sum of the radii 
of the ions one again has only elastic Rutherford scattering and 
at the most Coulomb excitation processes. 

(2) If b is comparable to the. sum of the radii of the ions, a 
grazing collision takes place and the ions can be elastically or 
inelastically scattered, or a few nucleons can be transferred from 
one to the other. Such collisions are collectively referred to as 
"peripheral collisions". 

(3) When b is reduced still further, the ions begin to interact 
very strongly, with a rapid loss of their energy and exchange of a 
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elastic scattering 
direct reactions 




elastic (Rutherford) scattering 
Coulomb excitation 


FIGURE 1 - 1 Schematic diagram showing heavy-ion trajectories at 
different Impact parameters corresponding to 
different reactions, (Ref, 18) 
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large number of nucleons between the target and the projectile. 
Because of the strong damping of energy from the relative motion 
these processes are. called deep-inelastic scattering. Such 
processes take place during the time intervals characteristic of 
the direct-reaction processes. In these reactions no fused 
systems are formed and the final product nuclei resemble the 
initial ones, but in highly excited states. 20 

(4) At still smaller values of the impact parameters, if b becomes 

equal to a critical value (b ) for fusion to occur, then after 

losing much of their energy in the relative motion the two nuclei 

form a joint system that rotates like a dumb-bell and further 

energy is lost in mutual excitations. At impact parameters below 

this value the two nuclei fuse together, completely losing their 

relative kinetic energy and their individual identities. The 

highly excited fused system then loses its energy by particle 

evaporation and/or by fission processes after a time interval 

which is much longer compared to the interaction times involved in 

21-23 

deep-inelastic scattering or direct reaction processes. 

At intermediate energies (of the order of 20-200 MeV/nucleon) 
the reaction processes are deep-inelastic collisons, incomplete 
fusion and fragmentation of the combined system into various 

n A 

nuclei. In near central collisions (b a* 0) at these energies, 

more violent reactions take place characterised by an almost 

complete disintegration of both the target and projectile nuclei. 

When the incident energy per nucleon is above the pion mass ^ 140 

3 24 25 

MeV, pion production may become increasingly appreciable. 
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In very high ultra-relativistic heavy-ion collisions (a few tens 
) the nuclear densities may rise so much that the 
nucleons overlap considerably, losing their identity and probably 
forming a Quark-Gluon-Plasma -experimental searches for which are 
still continuing .* 3 '^ 7 


1.3 THEORETICAL APPROACHES AND MOTIVATIONS FOR THE PRESENT WORK 


The de Broglie wavelength of a heavy-ion is given by 


X 


1/2 


2 M E 


lab 


( 1 . 1 ) 


The wavelength of a typical heavy-ion projectile such as ^°Ca at 


E lab = 115 MeV " an aner8y that is 

just 

sufficient 

to overcome 

the 

40 40 

fusion barrier for the Ca+ Ca 

reaction - is 

about 

0.42 

fm. 

whereas the experimental rms radius of 

^°Ca is 

about 

3.53 

fm. 

Thus heavy- ions with energies 

even 

of the 

order 

of a 

few 


MeV/nucleon have short de Broglie wavelengths compared to the size 
of the ions. Therefore, classical approximations are expected to 
be good at least for the macroscopic features of heavy-ion 
reactions, such as fusion and deep-inelastic scattering, which 
dominate at low incident energies and involve large scale transfer 
of energy from the relative motion to internal excitations of the 
colliding nuclei . 20 ' 23 Therefore, it has been common to use 

classical macroscopic approaches to study the collective motions 

, . . .... 20-23,26-36 

in heavy-ion collisions. 
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The macroscopic approaches are based on the assumptions of a 
continuous distribution of nuclear matter. In these calculations, 
therefore, one has to make a-priori choice of the relevant 
collective degrees of freedom describing the overall properties of 
the nuclei, such as separation between the ions, their shapes, 
neck elongation etc. After chosing the relevant collective 
degrees of freedom one has to invoke suitable mechanisms, such as 
friction, for the transfer of energy from the collective degrees 
to the remaining internal degrees of freedom describing the motion 
of nucleons. 

On the other hand, if one adopts a microscopic approach, 
which necessarily involves individual particle degrees of freedom, 
this energy transfer is determined automatically through the 
solutions of the time -dependent many-body equations; the relevant 
collective degrees of freedom are determined by the system itself. 
Exact solution of a quantum-mechanical many-body calculation 
would, of course, be ideal. But, even the best available 
approximate solution, such as the Time -Dependent Har tree -Foe k 
(TDHF) calculation 37-41 which is based on the mean-field 
assumption and a single Slater determinant wave-function, poses an 
enormous computational problem due to the non-linear integro- 
differential nature of the equations involved. Practical 
calculations require further assumptions to be made and several 
symmetries to be imposed on the system. Furthermore, as the 
incident energy is increased two-body collisions become 
increasingly important and the TDHF approximation becomes invalid. 
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It. has been felt, necessary to incorporate two-body collisions in 
TDHF calculations even at relatively low energies^ - such 
calculations pose even more computational difficulties. Therefore, 
classical microscopic approaches have also been considered which 
require orders of magnitude less computational effort compared to 
the TDHF calculations and are relatively clear in their 
interpretations . 

There are two types of classical microscopic calculations 

based on different approaches. One of the classical-microscopic 

approaches is based on classical kinematics and Monte-Carlo 
43”"47 

methods. In this approach, which is often called 

Intranuclear Cascade model , the nucleons are assumed to follow 
straight-line trajectories between scatterings. The result of 
each scattering is determined by Monte-Carlo methods according to 
the distance of closest approach of the two nucleons and the 
experimental values of the free N-N scattering cross sections. 
This approach necessarily requires the conditions of diluteness of 
the nuclear matter, such that the simultaneous interaction between 
more than two nucleons is explicitly excluded. The N-N 

interaction is replaced by the N-N scattering cross sections. 
Most of such calculations neglect binding energy effects and their 
application has been restricted to high energy collision 
calculations only, where two-body collisions become important and 
the binding energy can be neglected in comparison to the total 
energy. In similar Monte-Carlo framework, mean-field and Pauli 
blocking effects have been included recently through 
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the Vlasov-Uehling-Uhlenbeck approach .^ 6 66 

The second classical-microscopic approach is similar to the 

Molecular Dynamics method for the study of classical fluids and 

54-59 i . 

solids. In this classical -microscopic approach nucleons in 

the colliding nuclei are considered as classical, spinless, 

point-particles and the simultaneous interaction between all pairs 

, „ „ 60-77 

of nucleons is allowed through suitable two-body forces. 

The effective two-body potential is chosen either to give good 
free N-N scattering cross sections or to yield reasonable 
properties for the nuclear ground-states. In a collision 
calculation, trajectories of the constituent nucleons are computed 
on the basis of the classical non-relativistic equations of motion 
and the forces between all the pairs of nucleons. The nucleon 
trajectories thus computed can be used to study any desired 
macroscopic properties of the collision. 

While in this classical-microscopic (equations -of -motion) 
approach some of the limitations of the purely hydrodynamic 
approach 16 , such as, those of a local thermal equilibrium are 
avoided. At the same time hydrodynamics is automatically included 
as a limit. 60 ' 61 Further, the diluteness condition assumed in the 
Monte-Carlo type classical-microscopic approaches is also not 
required. Since, at higher energies (E lab >100 MeV/nucleon) the 
nucleon wavelengths (<1 fm) are comparable to, or smaller than the 
force range or the mean inter-particle separation, classical 
approximations are expected to be even better at such energies. 
Therefore, this approach has been extensively used in the study of 
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various aspects related to the formation and decay of the denser 

nuclear matter in high-energy (of the order of 100-800 

MeV/nucleon) heavy-ion collisions. Recently, this approach 

has also been used in the study of formation of composites in the 

69“7 1 

disassembly of "hot-drops" of nuclear matter. 


Some efforts have been made to use this approach at lower 
energies ( of the order of a few MeV/nucleon ) also. Since, all 
the translational degrees of freedom of the constituent nucleons 
are explicitly included in this approach, there is no need to 
separate the dynamics of the collective degrees of freedom from 
the intrinsic degrees, as in the case of classical-macroscopic 
approaches. Thus, there is no need to make any assumptions about 
the shape of the nuclear surfaces or the width of the neck formed 
between the two nuclei at contact. Further, there is no need to 
make any assumptions about the frictional forces, as all the 


dissipation and 

transfer 

of 

energy arises due 

•to 

the 

nucleon-nucleon 

collisions 

and 

are already built 

in 

the 


calculations . 


Ramamurthy and Kataria 72 used this approach for a 
qualitative study of heavy-ion collisions at lower energy and 
found that most of the essential features of heavy-ion collisions 
such as fusion and deep-inelastic scattering are adequately 


brought out in this approach. 
Ramamurthy and Kataria, Dixit et 


Following the calculations of 
al 73-75 and Mani 76 calculated 


fusion cross sections for various heavy-ion reactions in this 


approach. While these calculations qualitatively reproduce the 
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general features of the experimental fusion excitation functions 
and give good agreements on the higher energy side, fusion cross 
sections are overestimated on the lower energy side just above the 
Coulomb barrier . These deviations are attributed to the possible 
breakdown of the classical approximations at such low energies. 
On the other hand even a simple classical one-dimensional 
static-potential model^’ ^ ^ ’ ^reproduces fusion cross 

sections, precisely in this energy region - though such models 
fail at higher energies due to the increasing importance of 
dynamical effects which are neglected. Therefore, it is necessary 
to examine what is the real cause of the discrepancy between the 
experiments and the earlier classical-microscopic calculations. 

The details of the classical -microscopic approach and fusion 

cross section calculations are given in chapter 2 . Taking 

IB 1 fi 

. 0+ 0 reaction as an example, we re-examine Dixit et al's 

calculations in chapter 3 and try to explain the cause of the 
overestimation of fusion cross sections at lower energies in their 
calculations. The use of Lennard-Jones potential as an NN 
interaction by Dixit et al in heavy-ion collision calculations is 
also critically examined in chapter 3. 

Since a two-body potential and the ground-state nuclear 

properties are the central inputs to the classical-microscopic 

approach, it is desired to understand clearly how fusion cross 

sections depend on these inputs. While this approach allows for a 

systematic study of the repulsive -core of the effective NN 

66 67 

potential in high-energy heavy- ion collisions, * a systematic 
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study of fusion of heavy-ions at lower energies may provide 

information about its long-range attractive tail. Therefore, a 

systematic study of the sensitivity of fusion cross sections on 

the range of a soft-core gaussian NN potential is carried out in 

chapter 4. There we calculate fusion cross sections for 
40 40 

and Ca+ Ca reactions and compare them with the experimental 
data. In chapter 4 we also study the behaviour of heavy-ion 
collisions at higher energies. In particular, we examine the 
question of the existence of the 1-window for fusion which is 
predicted by the TDHF calculations, and present some calculations 
which demonstrate the "sideward-flow" or "side-splash" of nuclear 
matter at even higher energies. 

In chapter S we try to account for the finite size of 
nucleons, replacing the point nucleons of the classical- 
microscopic calculations by the "frozen" gaussian wave -packet s . 
Taking two different values of the rms radius of the nucleons we 
study the effect of the finite size on the classical NN potential, 
nuclear density distributions, compressibility coefficient and 
heavy-ion fusion cross sections. 

One of the difficult problems in the classical-microscopic 
approach is to explicitly incorporate the zero-point motion of 
nucleons in the nuclei and account for the "Pauli-repulsion" 
between the nucleons. Because of the classical nature of the 
approach and the quantum mechanical nature of the problem there 
does not seem to be any satisfactory solution to this problem 
within the framework of the classical-microscopic approach. 
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Therefore, in chapter 6 we finally describe a semi cl as si cal - 
microscopic approach with gaussian wave-packets which removes the 
“frozen-width" assumption on the wave-packets made in chapter 5. 
This approach explicitly incorporates the effects of the 
zero-point energy and the "Pauli -repuls ion" in the ground-state 
and in the dynamics of the wave-packets. Some simple calculations 
are also presented in chapter 6. 
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CHAPTER 2 


CALCULATIONAL DETAILS OF THE CLASSICAL-MICROSCOPIC APPROACH 


2.1 INTRODUCTION 


In classical-microscopic approach, nucleons are considered as 
point particles, which interact with each other via two-body 
forces. Appropriate charges are assigned to the particles. 
However, nucleons are assumed to be without any spin degree of 
freedom in the classical-microscopic approach. The motion of each 
nucleon is governed by the coupled Newton’s equations of motion : 


m 


, 2 + 

d r. 


dt 


-* 

-v 


i ( Z v ij } 


( 2 . 1 ) 


where, m is the mass of a nucleon which is assumed to have the 

same value for protons and neutrons; and v. . is the two-body 

i J 

potential which is chosen to yield reasonable ground-state 
properties of the nuclei. 

The different steps involved in the simulation of a heavy-ion 
reaction (A^+Ag) in this approach are 

i) Choosing a suitable two-body interaction potential, 
ii) Construction of the nuclei in their ground-state as initial 
configurations for the collison calculations. 
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iii) Specifying the initial conditions for collisions, 

iv) Finding the trajectories of all the A = A^ + Ag nucleons for 
the above initial conditions; and 

v) Analysis of the trajectory data to calculate reaction 

properties. 

In this chapter we describe in detail the above steps 
involved in the calculation of fusion cross sections. 
Calculational details are similar to those in the earlier 
calculations in Ref. 1-5. 

Procedures for determining the ground-state configurations of 
the nuclei are described in section 2.2. The necessary features 
of the two-body potentials in the present approach are discussed 
in section 2.3. Specification of the initial conditions for 
collisions and the calculations of the nucleon trajectories are 
described in section 2.4. Definition of fusion reactions and the 
details of fusion cross section calculations are given in section 


2.5. 
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2.2 CONSTRUCTION OF THE NUCLEI IN THE GROUND-STATES 

The individual nuclei used in the collision calculations are 
obtained by generating a random distribution of all the nucleon 
positions in a sphere of radius R, where R is approximately the 
radius of the corresponding nucleus. However, this random 
distribution of nucleon positions is generally very "hot'', i.e., 

energetically very unstable, since, in particular nucleons may be 
sufficiently close to each other for the short-range repulsion to 
be very large. Further, such configurations in general may not 
have density distributions that correspond to nuclear saturation 
densities. Therefore, to obtain a configuration of the nucleon 
positions corresponding to an equilibrium ground-state with the 
correct two-body correlations, this "hot" distribution may be 
"cooled" by one of the following two methods called the “DYNAMIC" 
method and the "STATIC" method, respectively. 

2.2.1 Algorithms 

Ca!> "DYNAMIC" Cooling Procedure s 

An obvious method of cooling a hot system is to put it in a 
"heat bath". Particles of the system collide with the "imaginary" 
particles of the bath and in the process lose their kinetic energy 
to these particles, while at the same time interact with other 
particles of the system. Practical steps involved in the 
calculations are as follows: 

(1) Velocities of all the particles in the system are initially 
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set equal to zero. 

(2) The “hot" distribution described above is then allowed to 

dynamically evolve following the classical equations of 

motion (eq. 2.1) with the chosen two-body potential v^ 

between the pairs of particles. 

-22 

(3) After a small time interval (usually^ 0.02 x 10 sec.) the 

particle velocities are again set equal to zero. This 

process removes the kinetic energy from the system and cools 
it; it also provides subsequent initial conditions for 

further evolution - the initial velocities are zero for the 
succeeding interval of the evolution. 

(4) At each predetermined time step various gross properties of 
the assembly of particles such as the total potential 
energy, and root mean square (r.m.s.) radius are monitored. 
The total potential energy is given by 

V = E E v (|r, - r . | ) (2.2) 

i j>i 1J 1 J 

where the sum is taken over all the pairs. The r.m.s radius 
is given by 

A -*■ ->2 

E (r. - R 0 ) 

2 1/2 r 1 1/2 

<x> L/c - 1^ J (2.3) 

where A is the total number of particles and r^and are the 
position vectors of the i lb particle and the centre of mass, 
respectively, with respect to the origin. 
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(5) The steps (1) through (4) are followed till all these 
macroscopic quantities settle down to their equilibrium 
values. Total potential energy V also acquires its minimum 
value at this stage. 

CbD ■■STATIC" method 

Another method to find a ground state configuration is to 
cyclically minimize the total energy of the configuration with 
respect to small displacements of single coordinates of different 
nucleons taken one at a time. We call this method "STATIC" because 
of the static displacements of nucleons. This method is similar 
to the static solution of the Hartree-Fock (HF) equations by the 
variational approach; in the present case wave functions being 
delta functions. Details of the algorithm are as follows: 

(1) The initial potential energy (U°) of an arbitrarily chosen i th 
particle is calculated 

°° = A Vij (|i ° - (2 - 4) 

and it is displaced from its initial position (r°) through a small 
step of size <5 along the x-axis, while keeping the positions of 
all the other particles fixed. 

Since all the other particle positions remain fixed, the 
total potential energy of interaction among these particles, viz 
(V - U°) also remains unchanged. 

(2) The potential energy of the i n particle in its new location 
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(r^) is calculated 

<--£ < 2 - 5 > 

and compared with the potential energy (U?) before it was 
displaced. If ^ 0°, then the above displacement is allowed. 

(3) But if > U°, then the particle is displaced along the 
x-axis by a step equal to 2<5 in the opposite direction. The 
potential energy of the l particle in this new location (r^ ) is 
calculated 

U i‘ = (2.6) 
and compared with U?. 

(4) The new position is allowed if U°, else the particle is 
restored to its original position (r?). Thus the particle is 
moved along a direction of decreasing potential energy. 

(5) The i tb particle is then moved along the y-axis from the final 
position of the above procedure and steps 1 through 4 are repeated 
for the y-coordinate of the i th particle; then similarly the steps 
1 through 4 are repeated for the z-coordinate of the i lh 
particle. 

(6) After the x,y, and z-coordinate searches have been made for 
the i lh particle the above procedure (steps 1 through 5) is 
repeated for the next particle. 

(7) This process is continued till all the 3A coordinate-searches 
fail in displacing any particle to a new location. At this stage 
the macroscopic properties of the cluster, like the total 
potential energy (eq. 2.2) and the rms radius (eq. 2.3) are 
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calculated and compared with the previous values. 

(8) The step length 6 is then reduced (usually by half the 
previous value) and steps 1 through 7 are repeated till the total 
potential energy converges to a minimum value. 

Since the “STATIC" method involves computations and 
comparisons of potentials only, it is found that when A is not 
very large, this method is much faster than the "DYNAMIC" method, 
which involves calculations of gradients of the potentials and 
numerical integration of eq. 2.1. A dynamical evolution of the 
final configuration obtained by the "STATIC" method confirms the 
dynamical stability of the configuration. Depending upon the 
final step size used in the "STATIC" method, only a small amount 
of total kinetic energy is found in the dynamical evolution which 
can be further removed by the use of the "DYNAMIC" cooling method, 
requiring only small lengths of time. 

It may be pointed out here that while for A=2, 3, and 4 one 

can find unique minimum- energy configurations, viz. an equilateral 

triangle, and a tetrahedron configurations for A= 3 and 4 

respectively, for A£5 a number of static isomeric configurations 

may exist. There is no practical method of determining whether a 

given minimum obtained numerically is in fact absolute or nearly 

6 

local in character. The only way out here, as is also done in HF 
calculations of the ground-state, is to take a large number of 
randomly generated different initial configurations and choose a 
final equilibrium configuration that has the minimum potential 
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energy among them. 

2. 2. 2 Frozen ground-states : 

The minimum potential energy configurations obtained through 
either of the methods in section 2.2.1 correspond to a frozen 
system at zero temperature when the particles are arranged 
approximately in a lattice with zero velocities. 

The frequently used technique of thermal excitation of the 
nuclei to simulate zero-point motion by arbitrarily assigning 
finite kinetic energies to the nucleons or even using a Fermi 

distribution of velocities results in an unavoidable premature 

7-12 13 

evaportion of nucleons. In Molecular Dynamics calculations, 

a standard practice to avoid the problem of evaporations is to 

impose periodic boundary conditions on the system so that a 

particle leaving the system from one end reappears from the 

opposite end. However, such boundary conditions cannot be imposed 

on the nuclei because of their finite size. 

While the. time scale for such evaporations is large compared 
to the collision times envolved in high-energy collisions, this 
time scale is clearly small compared to the time scales involved 
in low-energy collisions. In a low-energy collision process, 
therefore, premature evaporation of even a few nucleons would 
amount to extra dissipation of the relative kinetic energy of the 
colliding residual nuclei which seriously affects fusion cross 
sections. Therefore, no arbitrary assignment of kinetic energies 
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to the nucleons is made in the ground state nuclei. Thus the 
zero-point energy arising from the uncertainty principle and the 
Pauli exclusion principle is neglected explicitly in the 
calculations presented in this chapter. 


2.3 TWO-BODY POTENTIALS 


Due to the lack of the Fermi-gas zero-point motion as 
discussed in section 2.2.2, the internal kinetic energy in the 
ground state of a nucleus is zero and the binding energy is equal 
to the total potential energy. Therefore, the use of a free 
nucleon-nucleon (NN) potential would lead to particle clusters 
having very large binding energies per nucleon and also very large 
nucleon densities in their ground- states . Therefore, we use a 

purely phenomenological form for the NN potential and demand that 
it reproduce the main gross characteristics of the nuclei, such as 
the ground state binding energy and the rms radius etc. This is 

easily done by suitably choosing the constants of the effective NN 
potential such that the equilibrium macroscopic properties of the 
particle clusters in their ground state are close to those of the 
corresponding nuclei . 


The usual Coulomb potential between the charged particles 
( protons ) given by 


V C (r i - r J> = 


1.44 



(MeV) 


(2.7) 


is also added to the effective two-body potential. Further, the 
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differences in the effective strengths of the nuclear interaction 

between like and unlike particles is taken into account by taking 

the NN potential between the like particles to be about 20% weaker 

14 

’ than that between the unlike particles. Different forms of the 
NN potentials are considered in the following chapters. 

2.4 COLLISION CALCULATIONS 

• 2.4.1 Initial conditions 

We compute all the trajectories in the centre-of-mass 
reference frame. For practical considerations, it is assumed in 
all heavy-ion collison simulations that initially (at time t = 0) 

the two colliding clusters are separated by a finite distance 
(Rc^j) , such that the short range strong nuclear interaction 
between the two nuclei is negligible and the Coulomb interaction 
is also small, so that the dynamical deformation and the internal 
excitation of the ions is initially negligible and not yet 
introduced by the forces. A relatively large centre-of-mass 
separation (R^) of about 20 fm is found suitable for both the 
and ^®Ca+^Ca reactions considered in this thesis. 

For times t < 0, i.e. before the actual calculations are 
started, the two ions are considered as point particles with 
given charges, masses, centre of mass kinetic energy (E^) and 
impact paramater (b), and are assumed to be moving along the 
classical Coulomb trajectories they would have followed from far 
in the past (t « 0 ). Figure 2.1 depicts the specification of the 




calculations in the centre-of -mass frame. 
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initial conditions at t = 0. The positions (f^ , f^) and 

velocities (^, res P©ctive centre of masses of the two 

ions with a given centre of mass separation R^ at t = 0 are 

obtained according to the Coulomb trajectory calculations 
described in the Appendix A. 

At t=0 the two particle clusters in their ground state are 
placed at and respectively and all the nucleons are 

boosted by uniform velocities corresponding to the respective 
centre of mass velecities v^ and vv, (see Figure 2.1) . 

2.4.2 Trajectory Calculations 

With the above initial conditions the trajectories of all the 
nucleons are obtained by numerical integration of A coupled 
Newton’s equations of motion (eq. 2.1), as the two nuclei approach 
each other. 

For numerical solutions of the equations of motion (eq.2.1), 

each second order differential equation is split into two 

first-order differential equations and integrated using the 

15 

Fourth-order Runge-Kutta method. Small integration time steps 

(<5t) are required to ensure adequate conservation of linear 

momentum, angular momentum, and the total energy, while at the 

same time keeping to reasonable computing times. For 

-22 -22 

<5t=0. 02x10 -0.04x10 sec. the total energy is conserved to 

-22 

better than 1% and for 6t=0.1xl0 sec. it is about 8%. 
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2.5 FUSION CROSS SECTION CALCULATIONS 

2.5.1 Definition of Fusion 

Since fusion is commonly understood as a reaction that leads 
to a complete amalgamation of the two nuclei forming a compound 
system which is equilibrated in all its intrinsic degrees of 
freedom, it is necessary to understand clearly what is the 
experimental fusion cross section. 

When two heavy-ions "fuse" together, they usually form a 
compound nucleus or a complex which is similar to it, with some 
excitation energy and angular momentum. Consequently s it 
de-excites by emitting light particles and y-rays leading to a 
residual nucleus called “evaporation residue", and, if the fission 
barrier is small or reduced sufficiently by angular momentum, it 
fissions. In some cases the angular momentum of the fused nucleus 
is larger than the value for which the fission barrier of the 
compound-nucleus vanishes. In this case a compound nucleus is 
not formed but a kind of equilibrated two-centre system is formed 
which fissions into two fragments. Such a process is called 
"fast-fission" . 


The experimental fusion cross section (o^ ^ d e fi ne d 


as 


.16 


^fusion ^ER + ^fission 


( 2 . 8 ) 


where, is the evaporation residue cross section and it 
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corresponds to nuclei with a mass close to one of the compound 

nucleus, and a.. . is the fission or the "fast-fission" cross 

lission 

section, which corresponds to products which have a symmetric mass 
distribution around a mean value which is about half the 
compound-nucleus mass. 


In eq. 2.8, dominates when light compound nuclei are 

formed in the reactions, whereas <y M , , dominates when heavy 

fission 

compound nuclei are formed. When very heavy target-projectile 
combinations are involved, with the product Z^Zg larger than about 
2500 to 3000, fusion cross sections go to zero. At projectile 
energies larger than about 10 MeV/nucleon fusion is incomplete 
and fast light projectiles can be emitted in the very beginning of 
the reaction and the two remaining fragments can subsequently fuse 
together forming a compound nucleus with a mass smaller than the 
mass of the total initial system. 

In theoretical calculations, when the aim is to calculate 
only the fusion cross section (<yj ^ ) and the decomposition of 

^fusion ^- Ir k° different channels arising from the fused system 

is not required, it is sufficient to calculate the cross section 
for the formation of the compound or compound-like fused system 
and follow its evolution for times longer than the times 
sufficient for deep-inelastic scattering to occur. 


Therefore, we operationally define fusion as a process in 

which a combined nucleus-nucleus system evolves for a sufficiently 

-22 

long period of time (£30 x 10 sec), without undergoing prompt 
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fission. This definition is similar to that adopted in TDHF 
17 18 

calculations. ’ It is also similar to that adopted in the 

macroscopic calculations where fusion predictions refer to whether 

the system will be caught behind a barrier in the internuclear 

potential or not, without any consideration of the question of 

19 

the subsequent evolution of the long-lived trapped system. 

Figures 2. 2, 2. 3, and 2.4 illustrate the behaviour of a 
heavy-ion collision at different impact parameters for a given 
collision energy. The Figure 2.2 shows the C.M. trajectories of 
one of the ion in the ^0+^0 reaction at = 67 MeV and 

impact parameters equal to 0,4,6, and 8 fm. Figure 2.3 shows the 
time evolu^tion of the C.M. separation (R^) and Figure 2.4 shows 
the energy in the relative motion as a function of time for the 
same reactions. At b = 8 fm there is negligible loss of the 

energy from the relative motion and the scattering is elastic with 
the C.M. trajectories being unbound. At b = 6 fm, a large amount 

of energy from the relative motion is lost into the internal 

excitations. But the energy loss is not sufficient to keep the 

two nuclei together for long time and the two nuclei eventually 
reseparate leading to deep-inelastic scattering. At b = 4 fm, very 
large amount of energy from the relative motion is dissipated and 
the two nuclei fuse together. In this case the nuclei stick to 
each other and revolve around a common axis, R^ is bound with 
small oscillations. For a head-on collision (b = 0 fm), almost 

entire energy in the relative motion is lost, the two nuclei stop 


each other and fuse. 




PM ~ 66.8 MeV and different impact 
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FIGURE 2.3 

Relative distance (R^) between the centres-of-mass 
of two ions as a function of time for the 16 0+ 16 0 
reaction at = 66.8 MeV and different impact 

P a ^ameters . 



0 4 8 12 16 20 24 28 

t (ICT 22 sec.) 

FIGURE 2.4 

Time evolution of the kinetic energy in the relative 
motion for the reaction at Eq^= 

MeV and different impact parameters. 


66.8 
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2.5.2 Expressions for fusion cross sections 

20 

Heavy— ion fusion cross sections are given by 

"fusion = "* 2 * 0 ‘ 2i+1 > t < 2 - 9 > 

where is the transmission coefficient for the fth partial wave 
and is the reduced wave number. For simplicity we calculate 
fusion cross sections in the Sharp-Cutoff model. In this 
approximation it is assumed that for a given there is a 

critical value of the angular momentum (£ ) , above which no 

fusion takes place, viz. , 

1 for l < ! cr 
0 for t > f Gr 



where l is the highest partial wave to fuse. Thus, for a given 
cr 

En M fusion takes place for angular momenta between zero and t 

cr 

only and the fusion cross is given by 


^fusion 


n 


l 

cr 

z 

1=0 


(2f + 1) 


Summing the series, we get 

o' , = n K 

fusion 

^ n ’k 


2 

2 


{Z + 1 ) 2 
v cr 

Z (Z + 1) 
cr v cr 


( 2 . 10 ) 


Therefore, an expression for classical fusion cross sections 


is given by 
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o . = n b 

fusion cr 


( 2 . 11 ) 


where b is the critical impact parameter given by 

S-# JL 


b = X 4 t (* +1) 

cr cr' cr ; 


( 2 . 12 ) 


2.5.3 Initial relative random orientations and Statistical 
averagi ng 

It is not possible to find a spherically symmetric particle 
distribution for a cluster with a finite number of particles. 

p> 

Figure 2.5 shows projections of nucleon positions in a randomly 
oriented ground state clusters of 0 L It is clear from the Figure 
2.5 that the distribution of the nucleons in the cluster is not 

spherically symmetric. This fact can also be noted by finding the 

21 l^O 

principal moments of inertia of the cluster. For the/ cluster 

shown in the Figure 2.5 , principal moments of inertia are in the 

ratio 1.00 : 1.35 :1.49. Therefore, the cluster is not completely 

spherically symmetric. Therefore, it becomes necessary to find an 

average of b in eq. 2.11 over different trajectory calculations 
cr 

in which initially (at t=0) the two clusters placed at ^ and 

respectively^ are randomly oriented with respect to each other. 

The changes in the initial relative orientations do not alter the 

relative positions of the nucleons within the same clusters and 

the initial ground-state configurations remain invariant. 

Generation of pseudo-random numbers and a rigid rotation 

21 

transformation matrix are used to give random orientations to 
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each clusters. 

For 16 0+ 16 0 system (A^Ag =32), b cr in eq. 2.11 is averaged 

40 40 

over 12 Initial random orientations. Since Ca+ Ca system has 
(A 1+ A 2 = 80) more number of particles, effectively there is 
considerable averaging even for each orientation. Therefore, it 
requires proportionately' l« s n ™ ber of orientations for averaging 

b . Therefore, only six orientations are considered in fusion 

cr ’ 40 40 

cross section calculations for the Oaf Ca reaction in Chapter 

4. Calculations with twice as many orientations for both the 

reactions result in less than 2* differences in the respective 

average fusion cross sections. Each critical impact parameter 

,, . . , . . , „n accuracy of about 0.2 fm. 

(b ) is determined to an 

Ci 
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CHAPTER 3 

CLASSICAL-MICROSCOPIC CALCULATION OF FUSION CROSS SECTIONS FOR 
THE 16 0+ 16 0 REACTION WITH THE LENNARD- JONES FORM OF NN POTENTIAL 

3.1 INTRODUCTION 

In an application of the classical-microscopic approach to 
heavy-ion collision calculations at low projectile energies, 
Ramamurthy and Kataria* chose a Lennard- Jones form of the two-body 
potential for the interaction between the nucleons and studied 
16 0+ 16 0 reaction at E lab = 128 MeV. Their calculations showed 
that most of the essential features of heavy-ion collisions at low 
energies, such as fusion, deep-inelastic scattering and 
particle-exchanges are adequately brought out in the 
classical-microscopic approach. 

Following the calculations of Ramamurthy and Kataria, Dixit 
et al 2-4 also used a Lennard-Jones form for the NN potential given 
by 

v(x“i j ) = Ac [ - (^/r^) 6 ] (3.1) 

16 16 

and calculated fusion cross sections for 0+ 0 (Ref. 2); and 

16 O+ 40 Ca and 40 Ca+ 40 Ca (Ref. 3) reactions at various energies. 
In a similar approach Mani calculated fusion cross sections for 
16 0+ 16 0 and 12 C+ 14 N reactions, albeit with a different form of 
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the NN interaction potential (Ref. 5). These calculations 

qualitatively reproduce the shapes of the respective fusion 
excitation functions for all the above mentioned reactions and 
give good agreement on the higher energy side of the respective 
fusion excitation functions. However, it is observed that the 
calculated fusion cross sections for all the above reactions are 
systematically overestimated on the lower energy side of the 
fusion excitation functions, just above the respective Coulomb 
barriers . 

These deviations from the experimental data at lower energies 

have been attributed to the possible breakdown of the classical 

2 3 

approximations. ’ On the other hand, even a simple classical 
one-dimensional static-potential model reproduces fusion cross 
sections, precisely in this energy region -though such models fail 
at higher energies due to the increasing importance of the 
dynamical effects which are neglected. Therefore, it is necessary 
to examine the real cause of the discrepancy between the 
experiments and the earlier classical -microscopic calculations. 

In this chapter we re-examine Dixit et al ’ s calculations 
1 fi If? 

taking the 0+0 reaction as an example and try to explain the 
cause of the discrepancy in classical fusion cross sections at low 
energies. The use of Lennard-Jones potential as an NN interaction 
in heavy-ion collision calculations is also critically examined. 
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3.2 IMPORTANCE OF THE CORRECT INITIAL CONDITIONS 

It is noticed that in all the above mentioned calculations, 
initially at time t=0 for collisions the centres of mass of the 
colliding nuclei, placed at a relatively large separation of about 
20 fm are, however, assumed to be moving along unperturbed 
straight line trajectories ( see Figure 3 of Ref. 2 ) instead of 
the Coulomb trajectories followed from far in the past as 
described in section 2.4.1 and Appendix A. Initial velocities at 
t=0 correspond to a given initial collision energy (E^) and a 
given impact parameter (b). Thus, in all the above calculations 
the important effect of the repulsion due to the Coulomb 
interaction between the two ions is neglected for t < 0 at large 
separations . 


The Coulomb potential between two ions of charges and Z 
separated by a distance R is given by 


V C (R) 


1.44 Z ± Z 2 
R 


(MeV) 


(3.2) 


16 

The Coulomb potential between two 0 ions separated by 20 fm is 

about 4.6 MeV. This Coulomb potential, though small compared to 

2 

the highest energy (E CM = 67 MeV) for which Dixit et al calculated 

16 16 

fusion cross sections for the 0+ 0 reaction, can not be 

neglected at lower energies ( e.g. at E CM =16.7 MeV). 
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Figure 3.1 shows the centre -of -mass (C.M.) trajectories of 
one of the ^®0 ions in the ^®0+^0 reaction at £-,..=16. 7 MeV and 
b=5.6 fm, with 

(A) the correct initial conditions as described in section 2.4.1; 
and 

(A’ ) with the neglect of the Coulomb repulsion between the two 
ions for t<0, as in Ref. 2. 

The two cases, marked by A and A’ respectively, show remarkable 
differences in the C.M. trajectories. For a given energy and an 
impact parameter, while the correct specification of the initial 
conditions lead to scattering of the ions (A), the neglect of the 
Coulomb repulsion in the specification of the initial conditions 
leads to a case of fusion (A’). Thus, it is due to the neglect of 
the Coulomb repulsion for t < 0 that, in Ref. 2, the trajectories 
with even higher impact parameters could lead to the fusion of the 
colliding nuclei for a given collision energy, giving rise to the 
higher values of the critical impact parameters (b ). Since 

U JET 

fusion cross sections are given by eq. 2.11, the higher values of 

b result in very high values of the fusion cross sections for 
cr 

16 16 

the 0+0 reaction at low energies. This effect is even more 
pronounced in the reactions that involve heavier ions, as in 
Ref . 3 , 4 . 

Therefore, to illustrate the above effect, we re-calculate 

16 16 

fusion cross sections for the 0+0 reaction, choosing the same 
form of the two-body Lennard- Jones potential for the interaction 
between the nucleons and the same values of the parameters c and <y 




the initially straight line trajectories for the two energies 
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taken from Ref. 2. The calculated fusion cross sections for 
£(-,^=16.7, 37.6, and 66.8 MeV are shown in Figure 3.2 for both the 
cases. Fusion cross sections calculated with the correct initial 
conditions as described in section 2.4.1 are denoted by LI in the 
Figure 3.2 and the fusion cross sections calculated by neglecting 
the Coulomb interaction for t<0 , as in Ref . 2 , are shown denoted 
by LI’. It is clear from Figure 3.2 that the fusion cross 
section with the initially correct classical Coulomb trajectories 
(LI) is considerably reduced at a lower energy (£^=16.7 MeV) 
compared to that with the initially undeflected straight line 
trajectories (LI’). 

At higher energies the effect of the Coulomb repulsion at 
large ion-ion separations is considerably reduced and, as shown in 
the Figure 3.1, the ions initially move along nearly straight 
line trajectories for E^=66.8 MeV in both the cases (denoted by B 
and B’). Further, the difference between the trajectories in the 
two cases at this energy is also very small. Therefore, the 
assumption of an initially straight line trajectory in the earlier 
calculations is good only at higher incident energies and fusion 
cross sections in both the cases (LI and LI’) do not show any 
appreciable differences at £^=66 . 8 MeV (see Figure 3.2). 




1 fi 1 R 

Calculated fusion cross sections for the 0+ 0 reaction: 

□ - potential LI and initially Coulomb trajectories; O 
potential LI and initially straight line trajectories 
(LI’); + -potential L2, which reproduces ground-state 
properties of 16 0. Calculated points with different 
potentials and initial conditions are connected to guide 
the eye. Experimental data are shown by o (Ref. 6), 
X (Ref. 7), A ( Ref . 8 ) , # ( Ref . 9 ) , ■ (Ref. 10), and V (Ref. 11). 
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3.3 DEPENDENCE OF FUSION CROSS SECTIONS ON THE GROUND-STATE 
NUCLEAR PROPERTIES 

The parameters c and a of the Lennard- Jones potential (eq. 
3.1) from Ref . 2 are also given table 3.1 under the set called LI, 
along with the calculated binding energy and the rms radius of 
0. The rms radius of D 0 obtained in the present 

calculations is slightly smaller than the value obtained in the 
Ref. 2 for the same set of parameters (see table 3.1). Therefore, 
due to a slightly smaller size of the ions in the present 
calculations fusion cross sections with the initially straight 
line trajectories (LI’) are also slightly smaller compared to the 
corresponding fusion cross sections calculated by Dixit et al in 
Ref .2. 

The rms radius obtained in the present calculation, as well 

as in Ref. 2 are, however, about 20% smaller compared to the 

16 

experimental rms radius of 0. Consequently, we find that fusion 

16 16 

cross sections for 0+ 0 reaction with the potential LI and 

with the correct initial conditions are also underestimated 

6 — 1 1 

compared to the experimental data. Thus, fusion cross 

sections depend on the ground-state properties of the nuclei (see 

also Ref. 12 ). Therefore, we readjust the parameters c and a such 

16 

that the experimental binding energy and the rms radius of 0 are 
reproduced correctly. The new set of parameters (called L2) are 
given in Table 3.1. The potential L2 is plotted in Figure 3.3 
and compared with the potential LI. The corresponding fusion 
cross sections with the potential L2 are shown in the Figure 3 . 2 
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TABLE 3.1 


Parameters of the Lennard-Jones potential (eq. 3.1) 

16 

calculated binding energies and rms Radii of 0 for 
calculations . 


and the 
different 




Dixit et al 2 

LI 

L2 

Experiment 

« (MeV) 


2.20 

2.20 

2.72 


o (fm) 


1.80 

1.80 

2.36 


Binding Energy 

(MeV) 

1 

CD 

O 

O 

-96.8 

-127.6 

-127.6 

rms radius 

(fm) 

2.24 

2.12 

2.73 

2.73 


106 £ 8 i. 





Leonard -Jones Potential 
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(denoted by L2). It is seen from the Figure 3.2 that at all the 
energies considered, fusion cross sections with this potential 
(L2) have very large values compared to those with the potential 

£? 1 i 

LI, as well as to the experimental data. 


3.4 SUMMARY AND DISCUSSIONS 

To recapitulate, Dixit et al obtained an overestimation of 

IB IB 

fusion cross sections for the 0+0 reaction at lower energies 

and speculated that this may be due to the possible breakdown of 

the classical approximations at low energies. It is however, 

shown above that this overestimation of fusion cross sections is 

rather due to the neglect of the Coulomb repulsion in the initial 

stages and incorrect specification of the initial conditions. 

With the correct initial conditions, on the contrary, we get an 

underestimation of fusion cross sections when the same set of 

parameters for the Lennard- Jones potential are used. This 

1 B 

underestimation is due to the smaller size of the 0 nucleus 
generated by this set of parameters (LI). Adjustment of the 
potential parameters to get the correct ground-state properties, 
however, results in an even larger overestimation of the fusion 
cross sections; notably at all collision energies. 

Calculated fusion cross sections with the potential L2 are 


almost uniformly 

overestimated over 

a wide 

range 

of 

collision 

energies (Figure 

3.2). Since it is 

expected 

that 

the 

classical 


approximations become increasingly good as the projectile energy 
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is increased, it is difficult to make any definite conclusion 
about the extent of the invalidity of the classical approximations 
only on the basis of the above disagreement of the calculated 
fusion cross sections with the experiments. On the other hand, 
it is also to be noted that the Lennard-Jones potential (eq. 3.1) 
has a very long-range attractive tail which goes like 4 £ 
large r^j, falling-off much slowly compared to a more conventional 
gaussian or an exponential form of attractive tail used in Nuclear 
Physics studies (Ref. 13). This long-range attractive tail of the 
Lennard-Jones potential tends to increase the attractive part of 
the ion-ion potential and the radius of the fusion barrier (Rg), 
and decreases its height (Vg). Thus, according to the simple 

one-dimensional potential model (see Ref. 14 and chapter 4, eq. 4.2 
) , an increase in Rg and a decrease in Vg lead to an increase in 
the fusion cross section for a given collision energy. Therefore, 
the above over estimation of the fusion cross sections for a 
Lennard-Jones potential (L2), which correctly reproduces 
ground-state properties of a nucleus, is mainly due to its 
long-range attractive tail. This indicates that fusion cross 

sections are not only sensitive to the ground-state properties of 
the nuclei alone, but also to the details of the two-body NN 
potentials; in particular to the range of the NN interaction. 

A systematic study of the sensitivity of fusion cross 
sections on the long-range part of the interaction would require a 
two-body potential whose range can be adjusted without 
significantly affecting the ground-state properties of the nuclei. 
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Since the Lennard-Jones potential has only two parameters (e and 

o'), which essentially determine the position (r . ) and the depth 

mm 

(^in) of the minimum of the NN potential, it does not allow one 
to make such a systematic study. Hence, the Lennard-Jones 
potential which is extensively used in the molecular-dynamic 
simulations of classical fluids, does not seem to be a suitable 
two-body potential for use in the study of heavy-ion reactions in 
the present approach. A systematic study of the sensitivity of 
fusion cross sections on the long-range attractive tail of an NN 
potential is presented in chapter 4, choosing a different form of 
the two-body potential. 
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CHAPTER 4 


16 1 6 40 40 

CLASSICAL-MICROSCOPIC CALCULATIONS FOR 04 O AND Ca4 Ca 
REACTIONS WITH A SOFT-CORE GAUSSIAN FORM OF THE NN POTENTIAL 


4-1 INTRODUCTION 

The central input to the classical-microscopic approach which 

is described in chapter 2, is a two-body nucleon-nucleon (NN) 

potential. It is, therefore, desired to understand clearly how 

heavy-ion fusion cross sections depend on this input. While, this 

approach allows for a systematic study of the repulsive core of 

1 2 

the effective NN potential in high-energy heavy-ion collisions, ’ 
a systematic study of the fusion of heavy ions at low energies may 
provide information about its long range attractive tail. 

In chapter 3, it is shown that the use of the Lennard-Jones 
potential as an NN potential results in overestimation of fusion 
cross sections due to its very long-range attractive tail. The 
Lennard-Jones potential has only two parameters which essentially 
determine the position and the depth of the minimum of the 
potential and thus determine the binding energy and the rms radius 
of a nucleus. Lack of an adjustable range in the Lennard-Jones 
potential does not permit a systematic study of fusion cross 
sections on the detailed nature of an NN potential. Therefore, 
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for this systematic study a more sophisticated form of the NN 
potential is required, keeping the number of essential parameters 
to a minimum. 

In nuclear structure studies it has been common to choose an 
exponential or a Yukawa or a gaussian form for the attractive tail 
of the NN potential and a suitable form factor for the hard or 

Q 

soft repulsive core. For the present study we choose a soft-core 

4 

gaussian potential given by 

v ij <r ij ) = - y o (1 ‘ C/r iJ ) exp <- r « /r o 2) (4 - U 

where V Q , C, and r Q are the depth parameter, repulsive core radius 
and range parameter respectively. For a given value of rQ, the 
parameters Vq and C are determined by the nuclear ground state 
properties. The range of the NN potential can be varied by 
suitably adjusting the range parameter r Q . Keeping in mind the 
large overestimation of fusion cross sections with a long-range 
Lennard- Jones potential, a gaussian form of the tail is more 
appropriate as it goes to zero much faster than the Lennard- Jones, 
an exponential or the Yukawa forms of the potential. The shape of 
the repulsive core is not expected to play any significant role in 
the determination of the fusion excitation functions because of 
the relatively low projectile energies involved in the fusion of 
heavy- ions 




59 


A systematic study of the sensitivity of fusion cross 
sections on the range of the NN potential given by eq. 4.1 is 
carried out in section 4.2.1. For this systematic study, we 
calculate fusion cross sections for the 16 0+ l6 0 reaction with the 
different values of the range parameter r Q . Fusion cross sections 
for a heavier system, 40 Cat 40 Ca, are also calculated and compared 
with the experiments in section 4.2.2. Several features of the 
results and the effective NN potentials are discussed in section 

4.2.3. 

In this chapter we also study some behaviour of heavy-ion 
collisions at higher incident energies. In particular, in section 
4.3.1 we examine the question of the existence of the 4-window 
for fusion, which is Predicted by the TDHF calculations ^above^ a 
certain threshold energy. In section 4.3.2 we present a 40 Cat Ca 
collision calculation at E CM = 418 MeV/nucleon which shows the 
"sideward-flow" or the "side-splash" of the nuclear matter which 
has been observed experimentally. Finally, the conclusions are 

given in section 4 • 4 . 
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4.2 FUSION CROSS SECTION CALCULATIONS 


4.2.1 The 16 0+ 16 0 Reaction. 

Systematic Study Of The Dependence Of Fusion Cross Sections 
On The Range Of The NN Potential: 


To study the sensitivity of fusion cross sections on the 

range of the NN potential (eq. 4.1) we calculate fusion cross 

1 R 1 fi 

sections for the 1 0+ 0 reaction with different values of the 

range parameter r Q . Three different values of Tq (1.8, 1.4 and 

1.2 fm) are chosen. Parameters Vq and C are adjusted for each 
value of rg , such that the ground- state configurations in each 
case reproduce the experimental binding energy and rms radius for 
16 0 nucleus. For r Q = 1.8, 1.4 and 1.2 fm, the set of parameters 
are called potentials PI, P2 and P3, respectively, and are given 
in Table 4.1 alongwith the calculated ground-state properties for 
16 0 . 


1 6 

The ground-^tate configurations for 0 obtained by the 

methods described in chapter 2 and potentials PI, P2, and P3 all 

have the same value of the binding energy and the rms radius. 

Ratios of the principal moments of inertia for these 

configurations are given in Table 4.2. The nearly same ratios of 

the principal moments of inertia indicate that these 

configurations have similar mass distributions and shapes. Figure 

6 

4.1 shows the pair-distribution functions n i . e . the number 
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TABLE 4.1 


Parameters of the effective NN Potential (eq. 4.1) and the 

16 


calculated binding energy and rms 

states with different potentials. 

radius 

of 0 in 

the ground- 


PI 

P2 

P3 

Experiment 

r o (fm) 

1.8 

1.4 

1.2 


V q (MeV) 

142.0 

730.0 

3360.0 


C (fm) 

2.19 

2.28 

2.35 


Binding energy (MeV) 

rms radius (fm) 

-127.6 

2.73 

-127.6 

2.73 

-127,6 

2.73 

-127.6 

2.73 
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TABLE 4. 3 


Ratios of the Principal Moments of Inertia in the ground-states of 
16 

0 with different potentials. 


Potential 

: 1 

: X 2 : 

I 3 

PI 

1.00 

: 1.34 : 

1.44 

P2 

1.00 

: 1.35 : 

1.49 

P3 

1.00 

: 1.35 : 

1.49 
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0 with the potentials PI, P2 and P3 . 
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of pairs n(r^j) for a given pair separation distance (r^), for 
these configurations. Pair-distributions for the 16 0 

configuration with the potentials P2 and P3 are quite identical 
and that with the potential PI also has peaks similar to those for 
P2 and P3. 

Thus, all the above ground-state configurations have nearly 
identical structural properties. The two-body potentials PI, P2, 
and P3 are plotted in Figure 4.2. These potentials have nearly 
identical depths (-2.9 MeV), and positions of minima (2.6 fm) . 
These NN potentials mainly differ in their long-range attractive 
tails - potentials with larger values of Tq have longer tails. 
Nearly identical minima of the two-body potentials give rise to 
the identical properties of their corresponding ground- state 
configurations . 

From the above discussion it is clear that any differences in 

16 16 

the calculated fusion cross sections for the 0+0 reactions 

with the above potentials should then be due to the differences in 

the potentials in the tail region only. Fusion cross sections for 
16 16 

the 0 +0 reaction with the potentials PI and P2 are calculated 

for Eq. = 16.7, 37.6 and 66.8 MeV, and with the potential P3 for 

E cm = 13.5, 16.7, 20.2, 24.1, 28.3, 32.8, 37.6, 48.3, 66.8 and 

104.5 MeV (corresponding to v^ =0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 

-22 

1.5, 1.7, 2.0 and 2.5 fm per 10 sec. respectively). The 
calculated fusion cross sections are shown in Figure 4.3 as a 
function of E^. Comparison of the fusion cross sections 
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FIGURE 4.3 


Fusion cross section Vs. for the 16 0+ 16 0 reaction wil 
the potentials PI, P2 and P3 . For a given potential tl 
calculated points at different values of Eli are connect* 
with straight lines to guide the eye. -Experimental dal 
shown by o (Ref. 8), X (Ref. 9), A (Ref. 10), (Ref. 11] 
■ (Ref 12) and V (Ref. 13). 
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calculated with the different potentials clearly shows the 
sensitivity of the fusion cross sections on the range of the NN 
potential. Potentials with a longer range give higher values of 
fusion cross sections, and as the range r^ decreases, fusion 
cross sections decrease for all the energies. 

To understand the above results, we consider the static 

one-dimensional potential model, according to which fusion cross 

7 

sections are given by 

V B (*=0) 

'’fusion = " R I<'=°> [ 1 I 1 < 4 - 2 > 

CM 

where Rg(£ = 0) is the barrier radius and Vg(£ = 0) is the barrier 
height for a head-on collision (i.e. I or b = 0). This model 

gives a reasonably good description of fusion excitation functions 
on the lower energy side, just above the fusion barrier. The 
model, however, does not give good agreement on the higher energy 
side due to the increasing importance of the dynamical effects. 
Even so, the expression 4.2 is useful in understanding the 
sensistivity of fusion cross sections on the range r Q . 

Therefore, we find the barrier parameters Vg(b = 0) and Rg(b 
= 0) for the ^0+^0 reaction with the different potentials in the 
sudden-approximation. In the sudden-approximation, the ion-ion 
potential is calculated as a function of the centre-of-mass 
separation (R CM ), keeping the configurations of the two nuclei 
frozen. Figure 4.4 shows the static ion-ion potential for the 
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system with the potential P3 for one of the 
initial-relative-orientation at b=0. The barrier parameters, 
averaged over 12 initial orientations for the different potentials 
PI , P2 , and P3 , are given in Table 4.3. As the value of Tq 
decreases, the NN potential (eq. 4.1) becomes less attractive in 
the tail region (see Figure 4.2). Thus the attractive part in the 
ion-ion potential decreases with the decreasing value of Tq , and 
therefore, the height of the fusion barrier (Vg) increases and the 
radius (Rg) decreases. Therefore, according to eq. 4.2 fusion 
cross sections decrease with decreasing values of r^ for a given 
value of E CM' This qualitatively explains the sensitivity of 
fusion cross sections to the range of the NN potential. 

In Figure 4.5, we plot fusion cross sections with the 
potential P3, calculated with eq.4.2 and the barrier parameters 
given in Table 4.3. Comparision of this result with the results 
of the dynamical-calculations with the same potential, shows that 
fusion cross section agree reasonably well in the two cases in the 
lower energy region. At higher energies, however, the 
one-dimensional potential model fails to reproduce the observed 
trend for fusion cross sections as a function of This is 
because of the increasing importance of the dynamical effects at 
higher energies which are neglected in the derivation of eq. 4.2 
due to the sudden-approximation. It is remarkable that the 
dynamical calculation of fusion cross sections in the 
classical-microscopic approach satisfactorily reproduces the 
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TABLE 4.3 


16 16 

Average values of the static-barrier parameters for the 0+0 
system (b=0) with different potentials 


Potential 

r o 

(fm) 

V B 

(MeV) 

(fm) 

PI 

1.8 

9.29 

9.43 

P2 

1.4 

10.41 

8.47 

P3 

1.2 

10.80 

8.32 
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FIGURE 

4. 5 



Fusion 

cross 

section 

Vs. for 

the 

1 6 n 16 n 

0+ 0 reaction. 

The 

solid 

line 

is the 

prediction 

of 

the static 1-dimensional 


potential model (eq. 4.2) with, the barrier parameters for the 
reaction given in Table 4.3 for the potential P3. Dashed 
curve shows the prediction of the dynamical calculations. 
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experimentally observed trend for fusion cross sections. 

1 fi 16 

In the Figure 4.3 fusion cross sections for the 0+0 

reaction, calculated in the classical -microscopic approach, are 

0.10 

also compared with the experimental data. Fusion cross 

sections with r^ = 1.8 fm (PI) and 1.4 fm (P2) are overestimated 

compared to the experimental data; while fusion cross sections for 
16 16 

the 0+0 reaction with r^ = 1.2 fm (P3) show good agreement 
with the experiments . 

Ajn 40 

4.2.2 The Ca+ Ca Reaction 

16 16 

Since fusion cross sections for the 0+ 0 reaction, 

calculated with the potential P3 (^q = 1.2 fm), agree well with the 

experiments, we use the same potential for the study of a heavier 

40 40 

system. Fusion cross sections for the Ca+ Ca reaction are 
calculated using the potential P3. Calculated fusion cross 
sections for = 60, 82 and 94 MeV are plotted in Figure 4.6 and 

compared with the experimental determination of the fusion cross 

., 14-16 

section. 

40 

The ground-state properties of Ca generated with the 

potential P3 are given Table 4.4. The parameters of the potential 

P3 give about 5% less binding energy and a 12% larger value of the 
40 

rms radius of Ca compared to the experimental values. As shown 
in chapter 3, section 3.3, fusion cross sections depend also on 
the ground-state properties of the colliding partners. The bigger 
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TABLE 4.4 


Parameters of the potential (eq. 4.1) and the calculated binding 

40 

energy and rms radius of Ca in THE ground-states with different 
potentials 



P3 

P4 

Experiments 

r o (fm) 

1.2 

1.2 


V Q (MeV) 

3360.0 

1155.0 


C (fm) 

2.35 

2.07 


Binding energy (MeV) 

-326.0 

-342.1 

-342.1 

rms radius (fm) 

3.96 

3.53 

3.53 
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size of ions tends to enhance the probability for fusion due to 
the increased overlap between the ions. Therefore, in Figure 4.6 
we find that fusion cross sections for the 40 Ca+ 40 Ca reaction 
calculated with the potential P3 are also overestimated compared 
to the experimental data. 

The inability of this classical-microscopic model in 

simultaneously reproducing the correct ground-state properties for 

both the 16 0 and 4( ^Ca nuclei with the same set of the potential 

parameters may be due to the neglect of quantum mechanical effects 

which might be playing an important role in the structure of 

these doubly magic nuclei. Nevertheless, a deviation of about 12% 

in the rms radius from the experimental value does not seem to be 

too bad for a purely classical model with a very simple form of 

the two-body NN potential (eq. 4.1). It may be pointed out that a 

similar situation also arises for the ground-state properties in 

the Hartree-Fock calculations which use a simple density 

17 

independent interaction. 

The above situation then requires that the potential 
parameters be slightly adjusted to obtain correct ground-state 
properties of the colliding ions. This becomes necessary in view 
of the dependence of fusion cross sections on the ground-state 
nuclear properties. Therefore, the parameters V Q and C of the NN 


potential are readjusted, 

keeping r^ = 

1.2 

fm. 

such 

that the 

experimental ground-state 

binding energy 

and 

the 

rms 

radius of 

4 ^Ca are reproduced correctly. Parameters 

of 

this 

potential , 



FIGURE 4.6 


Fusion cross 

section 

Vs. E cm for 

the 

40 Ca + 40 C 

reaction 

with 

the 

potentials 

P3 

and 

P4 

Experimental 

data 

are 

shown by o (Ref 

. 14) 

, A 

(Ref 


15) and o (Ref . 16) . 
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called potential P4, are given in Table 4.4. 

The potential P4 is plotted in Figure 4.7 and compared with 

the potential P3. Since both the potentials reproduce the 

40 

binding energy of Ca to within a difference of 5%, they have 

nearly the same depth. But the potential P4 has a smaller radius 

of the minimum compared to that for the potential P3, because the 

parameters of the potential P4 have been adjusted to reproduce a 

40 

smaller rms radius of Ca compared to the rms radius generated by 
the potential P3. 

40 

Figure 4.8 shows the pair-distributions n(r^) for the Ca 

configurations with both the potentials P3 and P4. In both the 

cases the pair-distributions have similar peaks, but the peaks for 

the ^Ca configuration with the potential P4 are shifted towards 

smaller pair-separation distances ( r ^j) compared to the peaks of 

the pair-distribution for the configuration with the potential P3. 

Thus, both the configurations have similar internal structure, but 

the configuration with the potential P4 is compressed compared to 

1 6 

that with the potential P3. The pair-distribution of 0 with the 
potential P3 is also shown in the Figure 4.7. It is interesting 
to note that the first few peaks for the and ^Ca ground-state 
configurations with the same potential P3 have nearly identical 
peak-positions. 

40 40 

The calculated fusion cross sections for the Ca+ Ca 
reaction with the potential P4 for = 60, 71, 82, 94, 121, 151, 
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FIGURE 4.8 Comparison of the pair-distributions in the ground-state of 40 Ca 
with the potential P3 and P4. Pair-distribution of with the 
potential P3 is also shown. 


7 


221 and 328 MeV (corresponding to v CM =1.2, 1.3, 1.4, 1.5, ' 1.7 

-22 

1.9, 2.3 and 2.8 fm per 10 sec, respectively) are also show: 
in the Figure 4.6. Fusion cross sections with the potential P- 
are lower than those with the potential P3, and a reasonably gooc 
agreement with experimental fusion cross sections is obtained witi 
the potential P4 in the lower energy region. 

The data point of Doubre et al 14 for E CM =150 MeV is shown 

in the Figure 4.6 as a lower limit to the fusion cross section of 
40 40 

the Ca+ Ca reaction, because no measurements have been made of 
the cross sections for fission following the fusion. For > 

100 MeV, no other measurements have been reported. Therefore, a 

comparison of the calculated results with experiments is not 

l 

possible at higher energies. Calcinated fusion cross sections 

(Figure 4.6), however, do show the generally observed experimental 
1 8 

trend of decreasing fusion cross sections with increasing 

energy at high energies. Classical macroscopic calculations with 

1 9 20 

the proximity forces, and a classical mean-field calculation 

40 40 

also show similar behaviour for the Ca+ Ca reaction and 

considerably higher fusion cross sections at E^ = 150 MeV, 

compared to the data of Doubre et al. Thus, careful determination 

40 40 

of the experimental fusion cross sections for the Ca+ Ca 
reaction at higher energies, where fission cross section may 
contribute significantly to the fusion cross section in addition 
to the evaporation residue cross sections, is needed. 
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4.2.3 Disucsslons 

As discussed in the chapter 2, section 2.2.1, the nuclei i 
the present calculations do not have any internal motion in thei 
ground-state and the Fermi-gas motion of nucleons or th< 
zero-point energy arising from the uncertainty principle and th< 
Pauli exclusion principle is explicitly neglected. The explicil 
neglect of the zero-point energy in the nuclei is then directlj 
reflected in the effective NN potential used in order to obtain 
correct gross nuclear properties. 

3 

For normal nuclear densities (p= 0.165 nucleon/fm ) the 

average separation between the nucleons is about 2.4 fm. 
Therefore, the minima of the NN potentials are also located near 
this value (2.6 fm) , with a correspondingly large repulsive core 
(2.0 - 2.4 fm) (see Figure 4.2) which gives correct rms radius in 
the ground states of a frozen configuration of point nucleons. 
Since the kinetic energy in the ground state is equal to zero in 
the present approach, the total binding energy is equal to the 
potential energy. Therefore, only a shallow depth of about 2-3 
MeV is required for the effective NN potential to obtain correct 
binding energies. Thus, we find that the effective NN potentials 
in the model calculations have large repulsive cores and shallow 
depths . 

A large repulsive core of the NN potential creates an 
excluded volume effect around a nucleon and prevents close 
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collisions between a large number of nucleons of the two nuclei. 

The net result of this effect is then somewhat similar to that of 

the Pauli exclusion principle which does not allow the nucleons to 

21 

come close together. However, the saturation properties of 

nuclei are due to a much more complicated balance between the 
attractive nuclear forces and the repulsion created by the 
exclusion principle. The inability of the model to simultaneously 
reproduce the nuclear properties in different nuclei with the same 
potential parameters may also be due to the inclusion of the Pauli 
repulsion only in a very approximate way by choosing a larger 
repulsive core. 

The close collisions between a large number of nucleons of 
the two nuclei are prevented at lower energies due to the large 
repulsive core. Hence, we find that at energies just above the 
Coulomb barrier and at the critical impact parameter (b cr ) the 
two nuclei just stick to each other like mud balls and keep 
revolving around a common axis (see Figure 4.9), with a very few 
nucleons transferred between them. Therefore, in the present 
calculations, internal excitations and the resulting shape 
deformation is the predominant mechanism for dissipation of the 
energy from the relative motion. Some energy also goes in the 
rotational motion of the combined system at non-zero impact 
parameters . 

Apart from the dissipation of the relative kinetic energy, 
fusion of two heavy ions also depends upon a critical balance of 
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the&t "tractive and repulsive forces in the system. Since the 
effective NN potentials in the present calculations have large 
core radii, a short-range attractive tail is required to 


satisfactorily 

reproduce 

the 

experimental 

fusion 

cross 

sections . 

Thus, we find 

that 

an 

effective NN potential 

that 

correctly 

reproduces the 

gross 

nuclear 

properties 

and at 

the 

same time 


gives a reasonably good agreement with the experimental fusion 
cross sections is required to be quite a narrow potential with the 
range parameter r Q = 1.2 fm (potentials P3 and P4). 

The compressibility coefficient (K) for infinite nuclear 
matter is given by 


a 


K = k 


^E/A 


a k" 


k = k T 


(4.3) 


where kj, is the Fermi -momentum. For a finite nucleus, K is given 
by 


K = 


2 

R 


a 2 E/A 


a r‘ 



(4.4) 


We calculate the compressibility coefficient (K) for with the 
potentials PI, P2, and P3. We calculate the compressibility 
coefficients by small radial expansions and contractions of the 
corresponding nuclei around their ground-state configurations and 
calculating the binding energies of the compressed or expanded 



FIGURE 4.10 

rpv flUXl £ oi\ 

binding energy^ Vs. the r B s radius of the 
compressed or the expanded configurations of 16 0 
with the potentials PI, P2 


and P3. 
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4.3 HEAVY -I ON COLLISIONS AT HIGH ENERGIES 

4.3.1 ^-Window For Fusion 

An unexpected outcome of the TDHF theory is the prediction of 
the existence of an angular momentum window (-^-window) for fusion 
at high bombarding energies above a certain threshold value, above 
which fusion of heavy-ions takes place only between a lower limit 
(^) and an upper limit (^). The fall of the fusion cross 
sections at high energies in the TDHF calculations is then not due 
to a limitation in the angular momentum (■£ ) or the impact 

C r 

parameter (b cr> ) "but instead due to the opening of this low 
/-window- Several experiments have been performed to check 
this prediction, but so far no conclusive experimental evidence 
has been found for the existence of the /-window for fusion. 

In a pseudo-par^ticle simulation of the TDHF calculation, 32 
it is shown that the non-fusion of heavy-ions at high incident 
energies and low impact -parameters follows a coherent flow through 
motion of nucleons of one nucleus past the surface of the other 
nucleus. The TDHF calculations on some asymmetric systems, such 
as *®0+ 2 ^Mg (Ref. 26), also clearly show such transparency effects 
of passing of one nucleus through the other. This result is, 
therefore, mainly due to the long mean-f ree-path of nucleons in 
the nuclei and the suppression of the two-body collisions in the 
mean-field calculations which give rise to the increased 
transparency of the nuclei in near central collisions above a 
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threshold energy. 

On the other hand, two-body collisions are naturally built-in 
the classical-microscopic approach through the two-body 
potentials. The mean-free-path of nucleons is also very short in 
such calculations (zero in the frozen ground-states). Therefore, 
classical-microscopic approach does not show any such transparency 
effects. The TDHF caculations^ predict opening of an ^-window 

for the reaction at a threshold energy of E^ = 27 MeV. 

16 16 

The present classical-microscopic calculation of the 0+ 0 

reaction with the potential P3 does not show any ^-window even 
upto = 150 MeV. Figure 4.11 shows a head-on (b = 0) 

collision of 16 0+ 16 0 at = 150 MeV which leads to fusion of 
the two nuclei. For E CM 5 150 MeV the two nuclei fuse at all the 
impact parameters from zero to b . Figures 4.12(a-c) show the 
16 0+ 16 0 reaction at E CM = 200 MeV and different impact parameters. 
For E™, > 150 MeV, and at low impact parameters a number of 
nucleons evaporate prematurly from the highly excited combined 
system, leaving behind a residual nucleus which breaks up in a 
relatively short period of time. Such events lead to what is 
called "incomplete fusion". The number of nucleons evaporating 
prematurly decreases as the impact parameter is increased, but the 
two nuclei reseparate leading to deep-inelastic scattering (Figure 
4.12c). As the incident energy is further increased we find that 
more and more nucleons evaporate in the near central collisions. 
Therefore, we conclude that there is no ^-window for fusion in the 



£ 

<H 


8 ^ 

t> 


+ &>^+ 
++&£ o 

$£<***> 
Q . C=D 

t> t> 


+ % + «. 




> V* 

+ (s.“ 


<*> O 
+ -+> + Jr 

ftp* 


*v + 

+*++- 


4* 


■*>* 

+ j<> +* 

\* 


° ^ 


* * 

*J?N- 

?* 

<> 




4. S^4- 

*4 

V 

»&> 

> 


FIGURE 4.11 Same as in Figure 4.9 but for E^ = 150 MeV and 

iScic %x f^YA e. 30 fiY) x Ao -ffA . 



= 20 




O 



ID 

II 




90 






92 


classical microscopic approach. In an earlier calculation also 
it was concluded by Mani 5 that there is no ^-window in the 
classical microscopic approach. 

It is interesting to note that in some recent calculations 
34-37 

also, it has been observed that the inclusion of the two-body 

collision terms in the mean-field TDHF theory brings about 
additional damping of the relative motion. The ^-window for 
fusion disappears in such calculations with reasonable relaxation 

O Q 

times for the two-body collisions. Some other calculations have 
also reported the disappearance of the ^-window when the two-body 
collisions are incorporated into their calculations. 

4.3.2 Sideward-Flow Of Nuclear Matter 

Another interesting aspect of heavy-ion collisions at high 

39-42 

incident energies is the experimental observation of the 

collective flow of nuclear matter at finite angles in the 
transverse directions. The flow angle ( e p>) increases smoothly 
from 0° in the peripheral collisions to 90° in the central 
collisions (b=0). 

Classical-Microscopic calculations also show the 

“sideward-flow" or the "side-splash" of nuclear matter in the 

high-energy heavy-ion collisions. Figure 4.13 shows the time 

40 40 

evolution of a head-on collision (b=0) in the Ca+ Ca reaction 
with the potential P4 at = 418 MeV/nucleon. Figure 4.14 shows 
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six initial-orientations of the 
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the effect of averaging over different initial orientations by 
superposing the b=0 event for six different initial orientations 
on the same plot. The sideward-flow of nuclear matter in the 
direction nearly perpendicular to the beam axis is clearly seen in 
the Figures 4.13 and 4.14. At higher impact parameters the angle 
of the collective flow decreases from 90°. Figure 4.15 shows an 
event at b = 4.5 f m . 

Here, we are interested in only showing qualitatively the 

observation of the collective flow in the present calculations. 

The flow angle (©j.) can be computed from a sphericity analysis of 

43 

the kinetic energy flow tensor or from a transverse momentum 
44 

analysis 

The distribution of flow angles calculated by Molitories et 

al^ for the Nb+Nb reaction in a similar classical-microscopic 

model shows good agreement with the experiments. Flow angle 

45 

distributions in nuclear hydrodynamics calculations also agree 
well with the experimental data. However, due to the transparency 
effects as discussed in section 4.3.1, the mean field calculations 


such as TDHF 

do not show any collective 

flow 

in 

the 

transverse 

direction^®. 

Cascade calculations^ ’ 

also 

do 

not 

show this 


effect. A detailed study of the collective flow of nuclear matter 

can be useful for the determination of the equation-of-state of 
49 


nuclear matter 
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4.4 SUMMARY AND CONCLUSIONS 

In this chapter, we have presented some studies of 
40 40 

and Ca+ Ca reactions in the classical-microscopic approach. We 

choose a soft-core gaussian form of the NN potential, which allows 

us to suitably vary the range of the NN potential. With this 

potential a systematic study of the dependence of fusion cross 

sections on the range of the effective NN potential is carried out 
16 16 

for the ,0+ 0 reaction. This study shows that the calculated 

fusion cross sections are sensitive to the long-range attractive 

tail of the effective NN potential. In TDHF calculations also, it 

is found that the use of the different Skyrme interactions gives 

23-25 

rise to large differences in the fusion cross sections. In 

the present case, potentials with a longer range cause increased 
attraction between heavy-ions and give higher values of fusion 
cross sections. With a proper choice of the parameters of the 
effective NN potential, we find that the calculated fusion cross 
sections for the and ^®Ca+^Ca reactions agree reasonably 

well with the corresponding experimental fusion cross sections. 

A qualitative difference shown by the classical-microscopic 
calculations with the predictions of the TDHF theory is the 
absence of any ^-window for fusion, which is predicted by TDHF 
calculations. So far there is no experimental evidence also for 
the existence of any ^-window. Recent calculations which 
incorporate two-body collisions in the TDHF calculations also 
support the absence of any ^-window. This indicates that two-body 
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c °lii s i OIls are important in heavy-ion collisions even at incident 
energies of the order of a few MeV/nucleons . 
Classical-microscopic calculations also show the experimentally 
observed "sideward-flow" of nuclear matter in high energy 
collisions. The sideward-flow of nuclear matter in the heavy-ion 
collisions arises due to the strong transverse coupling between 
the nucleons with the short mean-free-paths - an effect that is 
completely absent in thejTDHF calculations due to the neglect of 
two-body collisions. 

Thus, we find that the classical-microscopic approach not 
only gives a good description of heavy-ion collisions at high 
energies, where the validity of the classical approximations may 
be obvious, but even at relatively low energies it provides an 
adequate description of the macroscopic features of heavy-ion 
reactions such as fusion. Classical -microscopic calculations do 
not require any additional assumptions such as those of frictional 
forces or shape parametrizations - assumptions which are . often 
required to be made in the classical -macroscopic approaches to 
account for the dissipation of the angular momentum and the 
kinetic energy from the relative motion. 

However, the lack of internal kinetic energy in the 
ground-state of the nuclei in the present calculations is 
reflected in the shape and strength of the effective two-body 
potentials for the point-like nucleons. Due to the explicit 
neglect of this zero-point energy we find that the effective NN 
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potentials have a large repulsive-core and a shallow depth in 
order to reproduce the ground-state properties of the nuclei. The 
range of the NN potential is also small, if it also reproduces the 
heavy— ion fusion cross sections. Due to the narrow width of the 
NN potential the nuclear compressibility coefficient has higher 
values. The problem of the finite size of nucleons and that of 
the zero-point energy are further addressed in the following 
chapters . 
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CHAPTER 5 

THE EFFECTS OF THE FINITE SIZE OF NUCLEONS 
CTROZEN" WAVE- PACKET CALCULATIONS) 

5.1 INTRODUCTION 

The calculations presented in the previous chapters show that 
the classical microscopic approach gives a good description of the 
macroscopic features of heavy-ion reactions over a wide range of 
energies. In this approach, however, nucleons are considered as 
point particles. With a point particle description of a nucleus 
we find that the density of nuclear matter in the nuclei is not 
continuous . It is also difficult to introduce the zero-point 
motion in the ground state of nuclei due to premature evaporation 
of nucleons from the nuclei. The NN potentials which reproduce 
correct ground state properties of a nucleus in Chapter 4 have 
large repulsive core radii and shallow depths; the range is also 
small if it also reproduces the experimental fusion cross sections 
for the heavy-ion reactions. These features of the NN potential 
result in a narrow width of the potential and very high values of 
the compressibility coefficients for the nuclei. 

On the other hand, we know from the experiments that nucleons 
do have finite size. Charge rms radius of protons is 
experimentally estimated to be about 0.8 - 0.9 fm.^ Thus^the size 
of a nucleon is not negligible compared to the average 
internucleon separation in the ground state which is about 2.4 fm. 
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Therefore, "to study what effects does the finite size of nucleons 
have on the classical NN potential, compressibility coefficients, 
nuclear density distributions, and heavy-ion fusion cross sections 
we consider nucleons as extended objects in the configuration 
space. While it may be relatively easier to consider nucleons as 
simply sharp spheres of given radii for this purpose, in the 
present calculations we, however, choose gaussian density 
distributions to describe a nucleon. An advantage of using the 
gaussian functions is its direct analogy with the minimum 
uncertainty wave packets which are also gaussians. In fact, a 
semi-classical wave-packet description is developed in the next 
chapter in which the zero-point motion of nucleons is explicitly 
included in the ground state and the width of the wave packets drs' 
also allowed to evolve. In the present chapter, in order to study 
the effects of the finite size and for the sake of simplicity we 
assume, however, that the width of the wave packets remains 
constant in time, that is, we assume that the wave packets are 
"frozen" . 

Further assumptions and calculational details are • given in 
section 5.2. In section 5.3 we choose two different size of the 
gaussian wave packets and study its effect on the NN potential, 
nuclear compressibility coefficient, nuclear density distribution 

-4 0 16 16 

for D 0 and fusion cross sections for the 0+ 0 reaction. 
Finally, summary and conclusions are given in section 5.4. 



5.2 CALCUL ATT ONAL DETAILS 


5.2.1 Description of a nucleon 


We consider each nucleon as an 
configuration space and describe it by 
gaussian density distribution given by 


extended object in the 
a spherically symmetric 


p(r) = 


1 


( 2 no 


2 ) 3/2 


exp (- 



- <r >) 2 



(5.1) 


The above distribution is normalised with 


' 3 (5.2) 

S p (r) dr = 1 

We chose the above gaussian density distribution (eq. 5.1) 

for a nucleon because of its analogy with the minimum uncertainty 

' , 2 
wave packets which are given by 

1 


- 


( 2 ,^) 3/4 


exp 


- (\'< r i>> 2 1 . 

r 5 <p i >-(V<V ) ] t6 - 3 > 

L Ao . 


where <r ± > is the centroid of the i th wave packet given by 


<r i > = S r i 0 i (r i )d"'r i 


= S r i P(r i )d J r i 


(5.4) 


The gaussian density distribution (eq. 5.1) is obtained by the 


relation 
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The rms radius (r^) of si nucleon described by eq. 5.1 or eq. 
5.2 is given by 

r N = [<r 2 > - <r> 2 ] 1/2 (5.6) 

where <r> is given by eq. 5.4 and <r 2 > is given by 

<r 2 > = /r 2 p(r)d 3 r (5.7) 

For the distribution given by (eq. 5.1) , r^j is given by 

r N = (3 o 2 ) 1/2 (5.8) 


5.2.2 Assumptions 


Since in the present chapter we are interested only in the 
study of the effects of the finite size of nucleons on the NN 
potential, we make the following simplifying assumptions: 

(1) Since a generalised wave packet in 3-dimension is given by 


P (x,y, z) = p(x) p (y) p(z) 


(5.9) 


with 


p(x) = 


- (x - <x > ) ‘ 


( 2^> 1/2 


exp 


2cr ‘ 


X 


and similarly for p(y) and p(z), we have already assumed in the 
equation 5.1 that the density distribution is spherically 
symmetric with 


a 


(5.10) 
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(2) We assume that all the nucleons in a nucleus have the same 
value of the rms radius (r^j), irrespective of their location in 
the nucleus and that it has the same value for both the protons 
and the neutrons. Therefore, we set 

= <^ = er (5.11) 

(3) We choose cr such that the rms radius ( r^) of each nucleon is 
approximately equal to the experimental rms charge radius of the 
protons. We choose a value of 0.8 fm for r^. We consider the 
case of a higher value of r^ also. 

(4) Finally, we assume that the width of each wave packet remains 
constant in time, i.e. 

a (t) = constant = a (5.12) 

The above assumption leads to the "frozen" wave packet 
approximation in which we picture nucleons as simply extended 
objects, which do not change their shape or size in time 
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5.2.3 Equations Of Motion 

According to the Ehrenfest theorem the equations of motion of 
the mean value- of the wave packets <r i > are formally identical to 
the classical equations of motion, 


m 


d r. 


dt‘ 


V£ ± v ij (? i - > 


( 5 . 13 ) 


except that the quantities which occur on both the sides of this 

3 

equation are replaced by their average (expectation) values. 
Therefore, for the motion of the centroids <i\>, we get 


m 


d 2 <r 1 > 


dt‘ 


= ~ V i <v ij (r i - > ] (5 ' 14) 


where ^ r i ~ r j ) > is the average potential between the i and 

the wave packets and is given by the double folding integral 

<v ij (r i' r j )> = "*i (r i> v ij (r i' r j> h (r i> d3r i d3r j 


( 5 . 15 ) 


or 

v ij < r r r j > d 3 r i d3r j •■■ (6 - 16) 

For a central potential v. .(r.-r .)', the effective potential <v. .> 

X J 1 J X J 

is a function of the separation between the two centroids only, 
viz | <r^> ~ <r j > l • 
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Because of the finite size of the wave packets, <v i j( r i _r j) > 
is in general not equal to ( | <r i >-<r j ) | > . The equality is 

maintained only in the special case when v, . is at the most 
quadratic in (r^r^) or, if o- is very small. The right hand side 

A. ‘U. 

of the equation 5.14 is the average force experienced by the 1 
wave packet due to all the other wave packets in the system. 

Numerical procedure for calculation of the trajectories of 
the centroids is similar to the description given in chapter 2 for 
point particles. 

5.2.4 Ground-State Properties Of Nuclei 

The ground-state of a cluster of "frozen” wave packets 
representing a nucleus is obtained by finding the ground state 
configuration of the centroids of the wave packets by the methods 
described in chapter 2 for the point particle calculations. 
Calculations of the ground-state properties is described below. 

(1) The nuclear-matter density p(r) at any given point r in space 
is given by 

^ 2 
p(r) = Z \<p, (r ) r 

i=l 1 
A 

= Z p.(r) (5.17) 

i=l 1 

where A is the mass number of a given nucleus and p^r) is given 
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by equation (5.1). 


A cluster of finite number of point nucleons, or wave packets 
does not have complete spherical symmetry. Therefore, for the 
purpose of calculating an average radial density distribution, we 
take an average of p(r) (eq. 5.1) evaluated at a large number of 
randomly chosen points in the volume of the nucleus at a 
particular radius (r). 


(2) The rms radius (R) of a nucleus is given by 

R = <r 2 > 1/2 


r J"p(r) r 2 d 3 r ,1/2 
«- /p(r) d 3 r J 


(5.18) 


where p(r) is the nucleon density given by eq.5.17. With the 
expression for p^ given by equation 5.1, we get 

A 2 2 

2 ( <r ± > + 3ap 

<r 2 > l/2 _ £ ill J 1/2 (5.19) 


With the second assumption (equation 5.11), we get 


<r > 


2 - 1/2 = [a, 2 + 1 y v 2 ] 


1/2 


and with equation 5.8, we get 


„ 2 . 1/2 
<r > 


A 1/2 

[4 + A ^ <r i >2 ] 


(5.20) 
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The second term In eq. 5.20 is the rms radius of the distribution 
of the centroids of the wave packets. The first term in eq. 5.20 
is the result of folding the wave packets on this distribution, 
which also corresponds to the size of a nucleus. 

(3) The binding energy E Q is given by 

A a- 1 A 

E q = Z <T > + Z Z <v (r.-r,)> (5.21) 

i=l i=l j>i 1J 1 J 

where, the first term is the total zero-point kinetic energy and 
the second term is the total potential energy. We have assumed 
here that in the ground state of a nucleus all the wave packets 
are stationary, i.e., 

<P i > = 0 

The zero-point kinetic energy <T^> associated with the wave 
packets is given by (see Chapter 6) 

, 2 V 2 
<P i > 3h 

<T . > = = ^ — (5.22) 

1 2m 2m (4crp 

Since we have assumed that all the o' s are constant in time, 
it follows from the eq. 5.22 that <T^> is also constant. 
Therefore, it contributes only a constant term to the binding 
energy (eq. 5.21) in the present case. Therefore, in the 

calculations that follow we first neglect this term in the 
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expression for the binding energy (eq. 5.21) and in the latter 
calculations we include this term also. 


(4) The nuclear compressibility coefficient (K) is given by 


* (E c /A) 


K = k. 


F 2 

*4 


For a finite nucleus it is given by 


(5.23a) 


P-f>. 


K = R‘ 


* (Eq/A) 


d R" 


C5.23>b) 


R=R 


g. s. 


4 

R is the rms radius in the ground state, 

g. s . 


where 
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5*2*5 Potentials and Forces 


For the present study, we choose the same form of the NN 
potential given by eq. (4.1 ) i.e. 



<r ij> 


-V Q a -C/r iJ ) exp (- r^j/rj?) (5.24) 


where, r, , = j r - r | 

J i j 

We also choose the Coulomb potential between the point protons 


given by 

v u (r ij> = 1 < MeV > 


(£. 25 ) 


and consider the potential between like pairs to be weaker than 
that between unlike pairs by about 20% as in the case of 
calculations in chapter 4. 


The double folding integrals (eq. 5.16) for the above 
potentials for the nucleon wave packets given by eq.5.1 are 
evaluated in the Appendix B. The result for NN potential (eq. 
5.24) with the assumption (2) gives 
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<v ±j tr i _r j )> 



r 


3 

o 


<r 2 + 4c 2 , 3 '' 2 
o 1 


C r 

o 

< r u > (^ 2 ) i/2 


erf 


x exp 


[ 



2 2 
(r“+4<0 


(- 

) 


« <r . ,> 

•e 


2a (r^+4 a 2 ) 1/2 - 


n 


(5.26) 


where 


<r 


ij 


> = i 


<r i> - <^>1 


(5.27) 


is the separation between the centroids of the i^* 1 and the j*'* 1 
wave packets, and erf(x) is the error integral given by 


erf (x) 



2 

-t at 
e 


(5.28) 


The result for the Coulomb potential (eq. 5.25) is given by 

1.44 


<r u > 


<v Ij ‘ r r r j» = 


<r ij> 


erf 


1 2a J 


(MeV) 


(5.29) 


In the limit of a ■* 0 the expressions (5.26) and (5.29) reduce to 
the expressions (5.24) and (5.25) respectively. 

Magnitude of the average force between the two wave packets 
<-V V i j> is easily found by differentiating <V i j> with respect to 
<r. •>. Differentiating eq. 5.26 for the NN potential with respect 

J 
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to <r^j>, we get 
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(5.30) 


Similarly for the Coulomb interaction from eq. 5.29 we get 


<r. .> 


-<r. -> 


c r 1 r S 1 1 r 'll 

< - Vv ij > = i 44 r 15—] — **p( — — j 

J *- <r ij > ^ ^ 4 <r ij > ^ 4c^ 


(MeV) 


(5.31) 
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5.3 RESULTS AND DISCUSSIONS 

We consider "two different values of o corresponding to two 
different values of r N (rms radius of a nucleon). 

5. 3. 1 The case of ■ 0. 8 f m 

First we choose a value of <y such that the rms radius of each 
nucleon (r^) is equal to the experimental rms charge radius of 
protons which is approximately equal to 0 . 8 fm. Since .all the 
wave packets are frozen ? the zero-point energy (eq. 5 . 22 ) is 
constant and since we are interested in studying the effect of the 
finite size of nucleons on the NN potentials, we first neglect 
this constant term in the calculation of the binding energy in eq. 
5.21 and in the determination of the potential parameters. Later 
in section 5.3.1(c) we include this term also in eq. 5.21 and 
study its effects on the NN potential and the fusion cross 
sections . 

In chapter 4 we found that the fusion cross sections for the 
16 0 + 16 0 and ^Ca+^Ca reactions, calculated by considering nucleons 
as point particles in the classical microscopic approach and with 
the NN potential given by eq. 5.24 agreed well with the 
experiments for a value of the range parameter (r o ) equal to 1.2 
fm (potential P3 and P4 in chapter 4). Therefore, for the present 
study also we choose the same form of the NN potential and 
r^ = 1.2 f m . 

We adjust the parameters V Q and C such that in the ground 
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1 6 

state of 0 the rms radius given by eq. 5.20 and the binding 
energy given by eq. 5.21 (with <T ± >=0) have values close to the 
corresponding experimental values for the ^0 nucleus. This set 
of parameters (called potential <W1>) is given in Table 5.1 along 
with the ground state properties of 16 0. 

Ca) Effects Of The Finite Size On : 

C i 3 Coulomb and NN potentials 

The Coulomb potential between two point like protons eq. 
(5.25) and its folding potential eq. (5.29) for r^ = 0.8 fm are 
plotted in Figure 5.1. Comparison of the two potentials shows that 
they differ significantly from each other only for <r i j > < 1.60fm. 

In Figure 5.2 the NN potential v? . (eq. 5.24) with the 

1 J 

parameters W1 and cr=0) is plotted and compared with the potential 

P3 of chapter 4. While the potential P3 for point particles is a 

narrow potential with shallow depth and a large repulsive -core 

N 

radius, we find that v. . (potential Wl) in the present calculation 

X J 

(derived for a nucleon size of 0.8 fm) is a much deeper potential 

with a smaller repulsive-core radius and a smaller radius of the 

* 

minimum. 

N 

The effective double folding NN potential <v^> (potential 

<W1>) is also plotted in Figure 5.2 and compared with the 

N 

potential P3 and v^(Wl). We notice that the potential <W1> also 
has a smaller repulsive-core radius compared to the potential P3. 
The tail of the potential <W1> is longer than the tail of the 
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TABLE 5.1 


Parameters of the effective 

1 6 

ground-state properties of 0 

NN potential 

and the 

calculated 

Potential 


<W1> 

<W2> 

P3 

P5 

Expt 

r N (fm) 


0.80 

1.40 

0.00 

0.00 


r Q (fm) 


1.20 

1.20 

1 . 20 

1.20 


C (fm) 


1.34 

0.99 

2.35 

2.23 


V Q (MeV) 


590.0 

495.0 

3360.0 

2090.0 


Binding energy 

(MeV) 

-127.6 

-127.6 

-127.6 

-127.6 

-127.6 

CE<r 1 > 2 ] 1/2 (fm) 


2.61 

2.34 

2.73 

2.61 


rms radius (fm) 


2.73 

2.73 

2.73 

2.61 

2.73 

K (MeV) 


440.0 

170.0 

995.0 

1015.0 




V (MeV) 
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FIGURE 5.1 

The effective Coulomb potential (eq. 5.29) for 
gaussian wave-packets with = 0.8 fm and 1.4 fm 
are compared with the Coulomb potential (eq. 5.25) 
between two point charges (r^=0 fm) 


/ 
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(A®W) A 


Comparison of the effective NN potential <W1> for gaussian wave-packets 
( r N - 0.8 £m) to a potential W1 with the same parameters but for point 
particles (r N = 0 fm). The potentials <W1> and W1 are compared with the 
potential P3 for point particles also. 
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potential P3 (see Figure 5.3 also). Thus, due to the overlap of 

the two wave packets, we find that <W1> is spread out in both the 

directions, i.e., it has a smaller core and a longer tail, 

compared to the potential P3 with the same range parameter r . 

o 

However, the depth of the potential P3 and <W1> are nearly 
equal because both the potentials reproduce approximately the same 
total potential' energy, which in the speical case of <T^>= 0 is 


also the binding energy. 

The minimum of 

the 

potential 

<W1> 

is 

located at <r. .> = 2.51 

X J 

fm which is less 

than 

the r . 

min 

of 

the 

potential P3 (2.60 fm) . 

The difference in 

the 

positions 

of 

the 


minima of the two potentials (about 0.1 fm) is nearly equal to the 

difference in the ms radius of 0 with the potential P3 and the 

rms radius of the distribution of centroids (corresponding to the 

1 6 

second term in eq. 5.20) in 0 with potential <W1>. 

Cii} Nuclear Compressibility Coefficient CIO 

We calculate the nuclear compressibility coefficient (K) 

1 6 

given by eq. 5.23b. The calculated value of K for O with the 

potential <W1> is about 440 MeV (see Table 5.1). Since <W1> is a 

wider potential compared to the potential P3, we find that the 

16 

compressibility coefficient for 0 with <W1> is much smaller than 

the value calculated with the potential P3. This value is close 

1*6 

to the value obtained for 0 in the HF and RPA calculations in 
Ref. 5. 


Citt » - Nu c lear D e ns i t y — D i l~strd-butlon- 
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CHID Nuclear Density Distribution 

Due "bo the finite size of nucleons, the frozen wave packet 

calculation shows remarkable differences with the point particle 

calculations in the nuclear density distributions. The average 

radial density distribution of nucleons p( r) for 16 0 with r„ = 0.8 

JM 

fm and the potential <W1> is shown in Figure 5.4. For a 
comparison with the point particle calculations, we show in Figure 
5.5a the distribution of point nucleons in with the potential 
P3. Whereas the point particle calculation shows a discontinuous 
density distribution in Figure 5.5a, the density distribution with 
the frozen wave packet calcuation in Figure 5.4 shows a continuous 
and a nearly smooth density distribution having a distinct surface 
region with a tail. 

16 

Figure 5.5b shows the distribution of centroids in 0 
with the potential <W1>. This distribution is similar to that 
shown for the point particle calculation with the potential P3 in 
Figure 5.5a, but is shifted towards the lower radii because of 
the smaller rms radius of the distribution of the centroids. The 
density distribution with the potential <W1> in the Figure 5.4 has 
a bump and wiggles on this bump in the region (r = 1.6 - 3.4 fm) 

where most of the centroids of the wave packets are located. The 
wiggles indicate that r^ = 0.8 fm is not large enough to generate 
a smooth density distribution in the present calculations. 


A comparison of p(r) 

in 

Figure 

5.4 

with 

the 

experimental 

6 

density distribution for 

16 0 

shows 

that 

p(r) 

in 

the present 



0.8 fm <W1> 
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calculation has a large value at very small radial distances. 
The reason for this higher density at smaller radii is the 
presence of a single nucleon in the equilibrium position near the 
centre of the nucleus in the present case (see Figure 5.5b) which 
contributes excessively to the average density. 

Civ3 Fusion Cross Sections 

We calculate the fusion cross section for the 
reaction with the potential <W1>. Due to the kind of 
approximations involved, the present calculations are similar to 
those of chapter 4 with the point particles replaced by the 
centroids of the frozen wave packets and the NN potential replaced 
by its double folding integral (eq. 5.26), but the rms radius 
given by eq. 5.20. 

16 16 

The calculated fusion cross sections for the 0+ 0 reaction 

with the potential <W1> at = 16.7, 32.8 and 66.8 MeV, averaged 
over 12 initial random orientations, are shown in Figure 5.6 and 
compared with the point particle calculations with the potential 
P3. The calculations of Chapter 4 showed that the fusion cross 
sections for heavy-ion reactions are sensitive to the tail of the 
NN potentials. Since the potential <W1> is wider and has a tail 
which is slightly longer than the tail of the potential P3 (Figure 
5.3) we find that fusion cross sections with the potential 
<W1> are also slightly higher than the fusion cross sections with 
the potential P3. 




Comparison of the calculated fusion cross sections 

1 C 1C 

for the 0+ 0 reaction with the potential P3 and 

potential <W1> with r^= 0.8 fm. 
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CfcO An Equivalent -Radi us Point Particle Calculation 

In some classical microscopic calculations'^ which treat 
nucleons as essentially point particles, the parameters of the NN 
potential are fitted to the rms radius given by eq. 5.20 to 
account for the finite size of nucleons. No other finite size 
effect on the NN potential is explicitly considered. Therefore, 
for the purpose of comparison with such a prescription, we again 
consider nucleons as point particles. We take the Coulomb 

potential given by eq. 5.25 and determine the parameters of the NN 
potential given by eq. 5.24 for r Q = 1.2 fm, such that the rms 
radius given by eq. 5.20 and the binding energy of 0 are close 
to the experimental values. Parameters of this potential (called 
potential P5) are given in Table 5.1. 

Potential P5 is plotted and compared with the potentials <W1> 

and potential P3, in Figure 5.3. Since the rms radius of the 

1 6 

distribution of centroids (second term in eq. 5.20) in O with 

1 6 

<W1> and the rms radius of 0 with the potential P5 are 
identical (2.61 fm) , we find that the two potentials have nearly 
identical positions of minima and nearly the same depths. 
However, the potentials P3 and P5 also have larger 

repulsive-cores compared to that for <W1>. The potential <W1> is 
also wider compared to both the potentials P3 and P5. 

1 6 

The calculated compressibility (K) for 0 with the potential 
P5 is about 1015 MeV which is much larger than the corresponding 
value for <W1>; and it is nearly equal to that for the potential 



P3 . The distribution of nucleons in ^^0 with P5 is shown in the 
Figure 5.5c which is similar to that with the potential <W1> in 
the Figure 5.5b and is shifted inwards compared to the 
distribution in the Figure 5.5a for P3. 

The calculated fusion cross section for the 16 0+ 16 0 reaction 

with the potential P5 at =16.7 MeV is compared with the other 

results at the same energy in Table 5.2. Since the rms radius of 

the distribution of nucleons (with r^ = 0 in eq.5.20) in 0 with 

1 R 

the potential P5 is less than the rms radius of 0 with the 
potential P3, we find that the value of the fusion cross section 
with P5 is also smaller compared to that with P3 and with <W1>. 

Thus, from the above discussion we observe that the main 
effect of the finite extension of the nucleons on the NN potential 
is the smaller core radius and an increased width of the 
potential . 

C c) Inclusion Of Zero-Point Energy 

For the r„ = 0.8 fm case , we now include the zero-point 
kinetic energy CT^ given by eq. 5.22, in the binding energy (eq. 
5.21). Since all the wave packets are frozen, all the o’s are 
constant. Therefore, <T^> is constant and inclusion of this 
constant term in eq. 5.21 represents merely a re-scaling of the 
total potential energy. Therefore, the rms radius of the 
distribution does not change and only V q in eq 5.26 is readjusted 
to give the correct binding energy. This new set of parameters 
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TABLE 5.2 

Fusion cross sections for the 16 0+ 16 0 reaction at E CM = 16.7 Me V 

with different potentials. 


Potentials 

°fus 

(mb) 

<W2+> 

1327 

<Wl+> 

1130 

<W2> 

1106 

<W1> 

838 

P3 

764 

P5 

680 
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(denoted by the potential <Wl+>) is given in Table 5.3. 

The potential <Wl+> is plotted and compared with <W1> in the 

Figure 5.7. The zero-point energy <T^> corresponding to r^ 

- 0.8 fm has a rather large value (72.8 MeV) . To compensate for 

this large positive term in eq. 5.21, the potential <Wl+> has a 

very large depth (V m ^ n = -26.6 MeV) compared to the depth of <W1> 

(V m in = -2.8 MeV), though the minima of the two potentials are 

located at the same separation <r^j>. At large values of <r^j>, 

the potential <Wl+> has a strongly attractive tail compared to the 

tail of <W1>. Since the effect of the repulsive (positive) 

zero-point kinetic energy term is included only in the binding 

energy and thus in the potential, this repulsive term does not 

appear in the dynamics in the present calculations. Therefore, 

due to the strong attractive long tail of <Wl+>, we find that the 

16 16 

calculated fusion cross section for the 0+ 0 reaction at = 

16.7 MeV with <Wl+> has a very large value compared to that with 
the potential <W1> (see Table 5.2). A proper inclusion of this 
repulsive <T^> term is discussed in chapter 6. 



TABLE 5.3 


The ground-state properties of 16 0 and the parameters of the 
effective NN potentials <W1> and <W2> modified to include 
zero-point kinetic energy in the binding energy. 


Potential 

<Wl + > 

<W2+> 

Expt 

* N (fm) 

0.80 

1.40 


r Q (fm) 

1.20 

1.20 


C (fm) 

1.34 

0.99 


V o (MeV) 

5520.0 

1840.0 


<T>/A (MeV) 

72.8 

23.8 


Binding energy (MeV) 

-127.6 

-127.6 

- 127.6 

[£<r. > 2 ] 1/2 (fm) 

2.61 

2.34 



rms radius 


<W2*>\<W1>\ I <W1+) 
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(A»W) A 
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5. 3. 2 The Case Of r ■ 1.4 f m 

Since the zero-point energy <T i > for r N =0.8 fm has a very 

large value (72.8 MeV), we now consider a value of r„ such that 

N 

o 

< 'T^> is of the correct order of magnitude. Therefore, we choose 
r N = 1.4 fm for which <T i > = 23.8 MeV. 

Firstly, we again neglect <T ± > in eq. 5.21 to study the 

effects of the increased size of the wave packets. With r„ = 

N 

1.4 f m and r = 1.2 f m we adjust the parameters V and C such that 
o o 

the calculated rms radius (eq. 5.20) and the binding energy (eq. 

1 6 

5.21) of the ground state of 0 are close to the experimental 
values. This set of parameters (called potential <W2>) is given 
in Table 5.1. 

Cal Effects Of The Fnite Size On : 

Cil Coulomb and NN Potentials 

The Coulomb potential (eq. 5.29) for r N = !' 4 fm is plotted 
in Figure 5.1' and compared with the Coulomb potential between 
point charges (eq. 5.25). This potential is also comapred with 
the Coulomb potential for r R = 0.8 fm. , in Figure 5.1. The 
Coulomb potential for r^ = 1.4 fm differs significantly from the 
point charge Coulomb potential for <r^> < 3.0 fm. 

The potential <W2> is plotted and compared with the potential 
<W1> in Figures 5.7 and 5.8. The potential <W2> has an even 
smaller radius of the core compared to the potential <lfl>, and due 
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to the increased size of the nucleons it is spread over much more 
than the potential <W1>. Thus a comparison of the potentials P3, 
<W1>, and <W2> (i.e., for r^ = 0.0, 0.8 and 1.4 fm respectively) 
in Figure 5.8 shows that, as r^ is increased, the corresponding 
potentials become wider with a smaller core radius and an extended 
tail. 

The a = 0 potential corresponding to the potential <W2> 
(called the potential W2) is shown in Figure 5.9 and compared with 
the potential W1 . The potential W2. is much deeper than the 
potential W1 and has a smaller radius of the core. 

Cil) Nuclear Compressibility Coefficient CIO 

The calculated value of the compressibilty coefficient (K) 

for 16 0 with the potential <W2> is about 170 MeV, which is less 

16 

than the corresponding value for 0 with the potential <W1>. 

Ciiil Nuclear Density Distribution 

16 

The average density distribution p(r) (eq. 5.17) for 0 with 
r N = 1.4 fm and potential <W2> is also shown in Figure 5.4. We 
notice from this Figure that due to the increased spreading of the 


wave packets 

this density 

distribution 

is 

smoother than 

the 

distribution 

with the potential 

<W1> 

and 

r N 

= 0.8 fm. 

The 

distribution 

of centroids of 

the 

wave 

packets 

in 16 0 with <W2> 

is 


shown in Figure 5.5d. From Figures 5.5(a-d) we note that with the 
higher value of o-(or r^), most of the centroids are located inside 
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(A*W) A 


Plot of the potentials W1 and W2 with the parameters of the respective 
potentials <W1 > and <W2> but for the case of o* = 0 
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a thin spherical shell on the surface. The tail of the density 
distribution, or -the spreading of the nuclear matter in this case 
is brought about by the spreading of the wave packets, whereas for 
the smaller value of a or for the point particle case the 
spreading at the surface is accomplished by the distribution of 
the particles itself. 


Civ) Fusion Cross Sections 

16 16 

The fusion cross section for the 0+ 0 reaction at = 

16.7 MeV with the potential <W2> is given in the Table 5.2. A 

comparison with the other results in Table 5 . 2 shows that it is 

higher than the fusion cross section with the potential <W1>, and 

it is very high compared to the calculated value with the 

potential P3 as well as compared to the experimental data. 

1 6 

Though, the rms radius of the distribution of the centroids in 0 
with the potential <W2> is 2.34 fm, which is less than that with 
the potential <W1> (2.61 fm), this high value of the fusion cross 
section is clearly due to the very long attractive tail of the 
potential <W2> (see Figure 5.8) which is much longer than the tail 
of the potentials <W1> and P3. 
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Cb!> Inclusion Of The Zero-Point Energy 

As in the case of = 0.8 fm, we include the zero-point 
energy <T^> in the binding energy (eq. 5.21) in the present case 
also. The parameters of the potential <W2+>, with the inclusion 
of <T^> in eq. 5.21, are given in Table 5.3. The potential <W2+> 
is compared with the potential <W2> and <Wl+> in Figure 5.7. The 
potential <W2+> is deeper than the potential <W2>, like the 
potential <Wl+> which is deeper than the potential <W1>, but since 
<T^> for = 1.4 fm is less than the zero-point energy for 
r^j = 0.8 fm (see Table 5.2), we find that <W2+> is less deep than 
the potential <Wl+> near the minimum. The two potentials, 
however, cross over again due to the increased spreading in <W2+> 
for r^ = 1.4 fm and thus the tail of the potential <W2+> is longer 
than the tail of <Wl+>. Therefore, the fusion cross section for 
the 16 0+ 16 0 reaction at E CM = 16.7 MeV with the potential <W2+> . 
given in Table 5.2 is higher than the corresponding fusion cross 
section with the potential <Wl+>. 
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5.4 SUMMARY AND CONCLUSIONS 

In the study presented above we replaced point nucleons by 
the frozen gaussian wave packets in order to account for the 
finite sise of nucleons in the nucleus. We studied the effect of 
the finite sise of nucleons on the NN potential, nuclear 
compressibility, nuclear density distribution, and fusion cross 
sections by choosing two different values of the width of the 
gaussian wave packets, corresponding to the rms radius of a 
nucleon equal to 0.8 fm and 1.4 fm respectively. 

In chapter 4 we found that with the point nucleon 
calculations the NN potential that reproduced the ground state 
properties of nuclei correctly had a very large repulsive core 
radius; and to reproduce the heavy-ion fusion cross sections also, 
it was necessary to choose a smaller value of the range parameter 
(r o ). This resulted in a short range, narrow potential with a 
very high compressibility coefficient (about 1000 MeV). In the 
present calculations, however, we find that with an increase in 
the size of the wave packets, the NN potential which reproduces 
the ground state properties of a nucleus has smaller repulsive 
core radius and the width also increases, which results in a 
smaller value of the compressibility coefficient (K). However, 
due to the long tail and the increased attraction between the 
colliding nuclei, the calculated fusion cross sections are 
enhanced . 

With point particles the nuclear density distribution is 
discontinuous with peaks occuring at the sites of the nucleons. 
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However, with the increasing width of the wave packets in the 

present calculations, we find that the density distribution 

becomes smoother. For r^ = 1.4 fm, we find that in the ground 
1 6 

state of- 0 nucleus most of the centroids of the wave packets are 
located inside a thin spherical shell on the surface which is 
reminiscent of the shell structure of the nuclei. We also find 
that r^ = 0.8 fm, which is equal to the experimental rms radius of 
protons, is rather small for the present calculations which 
results in some fluctuations in the density and a very high value 
of the zero-point kinetic energy (73 MeV) . 

The present calculations demonstrate that, as compared to a 

point particle description, the use of wave packets gives a better 

-n.ucAex\£ 

description of a ncul - cus and the NN potential. However, with the 
frozen width approximation for the wave packets in the present 
calculations, inclusion of the zero-point kinetic energy in the 
binding energy calculations makes the NN potential even deeper 
with an even longer attractive tail. Since this repulsive energy 
term is constant ( because of the frozen constant width ) , the 
increased attraction in the NN potential only enhances the 
heavy-ion fusion cross sections still further. In the next 
chapter we propose a semi -classical wave packet approach which 
attempts to remove some of the assumptions made in this chapter. 
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CHAPTER 6 

SEMICLASSICAL MICROSCOPIC APPROACH WITH GAUSSIAN WAVE PACKETS 


6. 1 INTRODUCTION 

It was noted in earlier chapters that due to the problem of 
premature evaporation of nucleons from the nuclei, it is difficult 
to explicitly introduce the zero-point kinetic energy of the 
motion of nucleons in nuclei in the classical microscopic 
approach. By replacing the "frozen" gaussian wave packets for the 
classical point nucleons, the zero-point energy was incorporated 
in the total energy in calculations presented in chapter 5 
sections 5.3.1(c) and 5.3.2 (b). However, since the wave 
packets were assumed to be "frozen" i.e. the width of the wave 
packets remained constant, the zero-point energy presented a 
constant term in the total energy which modified the effective NN 
potential. As it is a constant term, it does not appear in the 
dynamics of the "frozen" wave packets. Further, the widths of the 
wave packets were chosen with a given value which was assumed to 
be the same for all the nucleons in a finite nucleus. 

Therefore, in the present chapter we remove the above 
constraint of the same constant value for the widths of the wave 
packets and develop a semi-classical microscopoic approach which 
also explicitly incorporates the effects of the Fermi gas kinetic 
energy or "Pauli-repulsion" 1 in the dynamics of the gaussian 



144 


wave packets . The details of the approach to find the static 

ground- state configuration and the dynamic evolution of a 

many-body gaussian wave-packet system are given in section 6.2. 

As an example we study the dynamic evolution of a simple two 

wave-packet system for some given initial conditions in section 

6.3.1. In section 6.3.2 we find the ground state configuration of 
1 6 

0. Finally, the conclusions are given in section 6.4. 


6.2 CALCULATI ONAL DETAILS 

6.2.1 Nucleon Wave Packets 


We assume that each nucleon in a nucleus is described by the 

minimum uncertainty single-particle wave functions, which are 

2 

normalized 3-dimensional gaussian wave packets. 

A general 3-dimensional wave packet is given by 

*i (x ± , y i , z ± ) = <P ± (x i ) ^(y^ 0 i (z i ) (6.1) 


where , 


- (x. - <x. >) 




( 2 - 2 )1/4 

x i 


exp 


[ ^ + “ <V‘ 


( 6 . 2 ) 


and similarly for ^(y^ and ^(z^. 2 For the sake of simplicity, 
as in chapter 5, we choose spherically symmetric gaussian 
wave-packets with 

a - a - a -a (6.3) 

x y z 
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% 

and 

* i <?i> = exp [ 

. . . (6.4) 

The properties of the gaussian wave-packets (eq. 6.4). 
described in chapter 5, are again listed below. 


1. 

/ ^ ^i dx dy dz =1 


(6.5a) 

2. 

<r ± > = S <p * ( r ^ 

> ? i 

(r ± ) d 3 ri 

(6.5b) 

3. 

<P jL > = S ^ i *(r i 

-> 

) (-ibV i 

) <f> ± (r i ) d 3 r i 

(6.5c) 

4. 

? 2 
<rp = <r t > + 

Zo\ 


(6 . 5d) 

5. 

2 2 
<P i > = <P i > + 

h 2 

i JL 

4 <7.2 

X 


(6 . 5e) 

6. 

(Ar A ) 2 = <r 2 > - 

<v 2 - 

and 

( 6 . 5f ) 


(Ap ± ) 2 = <P 2 > - 

<P ± > 2 , 

with 



(Ar i ) 2 .(Ap i ) 2 = 

^ 2 
- fi 

4 

and 



(AXi) 2 ( A P Xt ) 2 

h 2 

4 

etc . 



- <r i >)‘ 




* * <V ( V <$ i > >] 


4 
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7. P ± (r) = |^ i (r )| 2 


(r i - 


<r i >)‘ 


(2na 2 ) 3/2 


exp 


(6.5g) 


2cr‘ 


Many Particle Nuclear Wave— Function and Nuclear 

Hamiltonian 


The many-body nuclear Hamiltonian is given by 


-h‘ 


H = Z 
i 


2m 


7 i + 


z 

i 


j>i 


v ij (r i 


r j ) 


( 6 . 6 ) 


Since nucleons are fermions, the A-particle nuclear 
wave-function is antisymmetric with respect to the exchange of 
particles. However, in the present approach we assume that the 
A-particle wave function y is the product wave function of the A 
single-particle wave functions given by eq. 6.4, i.e. 

A 

V - n (6.7) 

i=l 1 


The important effect of the exclusion principle is included 
in the present approach as described later. However, at this 
stage we assume the Hartree-ansatz (eq. 6.7) and find the 
expectation value of the Hamiltonian (eq. 6.6), which is given by 
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-h* 2 

<H> = £ < V > + E E < v 

i 2m i j> i 


ij (r i 


r j» 


( 6 . 8 ) 


where 


<v. . (r , -r .> 
U i J 


^ i (r i )<A. ( rj ) v i j(^ i -r j )^ i (r i )<^ j (r j ) d 3 r i d 3 r.j 


(6.9) 


>2. 


and <P^> is given by eq. 6.5e. Therefore, we get 



<P. > 2 


. 2 


r r 

1 

3 

h 

J_ i 

Li 

*r 

2m 

4 

2m 



] + Z Z 
-* i j> i 


<v iy 


(6.10) 


The first term on the right hand side of eq. 6.10 is the total 
kinetic energy of motion of the centroids of the wave packets. 
The second term gives the quantum mechanical part of the total 
kinetic energy, and the third term is the total potential energy. 
The nuclear density p(r) is given by 

... ^2 


- (r - <r i >)‘ 


P(r) ~ ? (9na 2 ) Z/2 ^ 

l ( 2 no ^ ) 


( 6 . 11 ) 
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6.2.3 Determination of a 


The expression 6.10 is the same expression as that used in 
chapter 5 to calculate the total energy. However, in chapter 5 it 
was assumed that all the ^ remain constant and moreover specific 
values of <y ^ were chosen which were the same for all the nucleons. 
At this stage in the present approach we now remove these 
assumptions for o in eq. 6.10 and determine the width o ^ of each 
wave packet from the following considerations. 


The second term on the right hand side of the expression for 
<H> (eq. 6.10) is the quantum mechanical part of the total energy. 
In the ground state of a nucleus we assume that <P^> = 0» that is, 
the centroids of the wave packets are stationary. Therefore, the 
total kinetic energy in the ground state is equal to the quantum 
mechanical zero-point energy, 

T _ V 5 - — < 6 - 12 > 

8M ' 1 4 2m 4 


The average kinetic energy per nucleon (<T>/A) in the ground 
state of a degenerate Fermi-gas is given by 


<T>/A = 


3 


5 2m 


3n 2 2/3 p 2/3 

2 J 


(6.13) 


where p is the constant density of the nuclear matter. To account 
for the inhomogeneity in finite systems, usually terms containing 
gradients of the density p (WeizsScker term) 3 are added to the 
kinetic energy density. However, for simplification we neglect 
such terms and calculate the kinetic energy of a nucleon in the 
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local density approximation by substituting p in eq. 6.13 by the 
actual density p(r=<r>) at every nucleon site. Therefore, the 
total Fermi-gas kinetic energy (Tp) is given by 


3tt‘ 


T = Z - 
F i 5 


-1 2/3 p 2/3 


2m 


J 


(r=<r i >) 


(6.14) 


In the ground state the total kinetic energy calculated from 
eq. 6.14 must, therefore be equal to the total kinetic energy 
calculated from eq. 6.12. In excited states, we have <P^> * 0 and 
the total kinetic energy is given by the first two terms in <H> 
(eq. 6.10) while the average Fermi kinetic energy in eq. 6.14 
requires corrections for non-zero temperatures in the excited 
states. However, as far as the quantum mechanical part of the 
total energy in the excited state is concerned we assume that it 
can still be obtained from either eq. 6.12 or, eq.6.14, i.e. we 
assume that = Tj,. 


By equating the right hand sides of the two equations (eq. 
6.12 and 6.14) we get, for each nucleon 

3 h 2 1 3 ^ r 'j 2/3 ^2/3 

4 2m <y 2 5 2m *■ 2 •* ( r-<r ± » 

or 


p(r=<r i >) 






3/2 


(6.15) 


where p(r = <r ± >) from eq. 6.11 is given by 
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p(r=<r ± >) = 


j (2no^) 3/2 


exp 


- (< ri > - <rj>)‘ 

2^? 


(6.16) 


Separating the j = i term in eq. 6.16, and equating eq. 6.15 and 
6.16 for p(r=<r i >), we get, 


1 


«»*) 3/2 


+ E 
j*i 


(2 no 


1 


2.3/2 

y 


exp 


- ( <r . > - <r,>) 

r - l — ,i 


2& 


j 


J 


, 1 ^ 3/2 

k A } 




(6.17) 


Denoting Rj, Rg, R A for <x ± >, <r g > , . . . . <r A > respectively, for 

simplicity,, and solving for cr ^ in eq. 6.17, we get an expression 
for c given by 


a 


i 




exp 


( 



(6.18) 


where 



(2 rr) 


3/2 



1/3 


(6.19) 


For i = 1,2,.... A, eq. 6.18 forms a set of A non-linear coupled 
equations for which can be solved by a self-consistent 
iterative procedure. Starting with any given initial values of cr ± 
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the solution of eq. 6.18 gives a unique set of values for the 
width of the wave packets, o^’s do not depend on the choice of 
the interaction potential but are entirely determined by the 
instantaneous state of a given conf iguration {R^t)} . Therefore, 
we find 

a i = (Vt). R 2 (t), R A (t))s o' i ({R i (t)}) (6.20) 

The <y^’s do not have explicit time dependence. The implict time 
dependence, however, does appear in cr'^ s due to the time 
dependence of the configuration {R^(t)} which is described in 
section 6.2.6. Though we have used the product wave function 
ansatz for the nuclear wave function (eq. 6.7) which explicitly 
neglects all the correlations, we have built-in some correlations 
through the above prescription and the dependence of the width 
(o^) of every wave packet on the entire configuration {^(t)}. 

Figure 6 . 1 shows some examples of configurations for the 

simple cases of A=2, 3 and 4. The values of the width a ^ are 

determined from eq. 6.18 for the 3-dimensional gaussian wave 

packets (eq. 6.4) which are centred at the points shown in the 

figure. Figure 6.1 demonstrates the configuration dependence of 

o^’s. The A = 2 case is further discussed in section 6.3.1 and a 

1 6 

ground state configuration of an A = 16 ( 0) nucleus is obtained 


in section 6.3.2. 




FIGURE 6-1 


Some configurations of the 3-dimensional gaussian wave-packets in the 
simple cases of A = 2, 3, and 4. The points shown are the sites where 
the centroids of the wave-packets are located and the numbers are the 
width (o') of the corresponding wave-packets in fm. * 7 the length a is equal 
to 1.0 fm. 
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6.2.4 Two-Body Potentials 

For the present study we choose the same form of the two-body 
potentials which are used in chapters 4. The NN-potential is 
given by 

= -V 0 (1 - C/r^) exp (- r^/r* ) (6.21) 

and the usual Coulomb interaction between protons is given by 

v ij (r ij } = 1 • 44/r i j (MeV) (6.22) 

The expectation value <v^> (eq. 6.9) of the above potentials are 
derived in the Appendix B. The result for the NN potential (eq. 
6.21) gives 



gives 
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(r ij }> 


1.44 


R 


erf 


ij 


R ij 

f (2c? 2 + 2^j) 1/2 i 


( ffleV ) 


(6.26) 


If o , i = Cj = <y then <v ij( r ij) > ( e< l- 6.23) and <v ij( r ij) > (©<1 6.26) 
reduce to the corresponding expressions ( 5 . 26 and 5.29) in chapter 
5. 


6.2.5 Ground State Solution Of Many-Body Systems 

In the ground state we have <P^> =0 as discussed earlier. 
Therefore, the expectation value of the Hamiltonian <H> or the 
total energy is given by 

3 h 2 1 

<H> = Z - — . — - + Z Z <v > (6.27) 

i 4 2m i j>i J 

To find the ground state configuration and the ground state 
binding energy we use the variational principle. We choose the 
trial single-particle wave functions {< i > by picking up a randomly 
generated configuration {R^} of the position of the centroids of 
the gaussian wave packets. To start with, we also choose {o^} 
randomly. However, since {c^} are uniquely determined for a 
given configuration {R^} by solving eq. 6.18, we do not treat {&^} 
as variational parameters in this approach. We then cyclically 
minimize the total energy of the system given by <H> in eq. 6.27 
with respect to small variations «5R i of the coordinates of the 
centroids of individual wave packets, taking one at a time and 
keeping the other centroids fixed. This method is similar to the 
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"STATIC" method described in chapter 2 (Section 2.2. (a)). 


However, an important difference arises between the method 
described here and the "STATIC" method described in chapter 2 
(section 2.2(a). In chapter 2, it was sufficient to calculate and 
minimize the total potential energy of the individual nucleons in 
each variation, since, due to the fixed positions of all the other 
nucleons the interaction energy between the fixed nucleons does 
not change. However, it should be noticed that in the present 
method, though we keep all the other {R^} fixed for each 
variation, due to the dependence of cr on the positions of all the 
centroids {R^}, a single variation of' the centroid coordinate of a 
given wave packet changes the values of o' for all the wave 
packets. Therefore, the interaction energy between the fixed wave 
packets also changes due to the dependence of each on a ^ and 
o j (eq. 6.23, 6.26). Thus, it is necessary in this case to 
minimise the total energy (eq.6.27) instead of minimizing only the 
total energy of each wave packet given by 


U 

i 


3 h 

4 2m 


a 


j* i 


< \i > 


Except for this departure the basic structure of the algorithm, 
however, remains much the same as that of the method "STATIC" in 
section 2.2(a). 


The ground state binding energy is calculated from eq. 6.27 
with the ground state wave functions, and the rms radius is given 
by 
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, 2 . 1/2 
<r > = 


I ( <i‘ i > + 3o^) 


1/2 


(6.28) 


6.2.6 Dynamic Evolution of a Many-Particle Gaussian Wave Packet 
System. 

If H c z.( r i» r 2 r A , P 1 ,P 2 > • • • -P A ) is the classical 

Hamiltonian of a many-body system then in the Hamiltonian 
Formulation of classical mechanics the equations of motion of 

5 

and p i are given by 

^1 * H cl (r l, r 2---- r A- f l, p 2'-V (6.29a) 

dt d 


and 

dPi d H cl (r 15 r 2 . ...r A , Pl> P 2 --.P A ) 


dt d r^ 


(6.29b) 


We assume that the gaussian wave packets remain gaussian at 
all times, though this is strictly true only for the potentials 
with upto quadratic terms only. With this assumption the dynamic 
evolution of the wave packets is obtained by the equations of 
motion for the centroids of the wave packets {R i } while the widths 
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{er^} are determined in the present approach by eq. 6.18 from the 
instantaneous configuration of the centroid positions (R^(t)}. 

We adopt a semi-classical picture in which we replace the 
classical variables {r^, P^} and the classical Hamiltonian 

H ci {r^, P^} by their respective expectation values, viz 
{R^ = <r^>, = <p^>} and <H> (eq. 6.10). For the two-body 

potentials which are chosen in section 6.2.4, it is clear from 
the equations 6.23 and 6.26 and from eq. 6.20 that <H> given by 
eq. 6.10 is a function of {R^} and{P^} only, that is 
<H> = S*f({R^,P i >) . Therefore, we also replace #/r^ and 9/d-p^ in 
equation 6.29 by tf/dR^and #/#P^ respectively. Thus, we get the 
equations of motion for the mean values {R^} and {P^} of the wave 
packets from the eq.6.29; 


dR^ dar({R i ,p.}) 

dt ' 3P k 


(6.30a) 


dP k «af({R i ,P i }) 
dt 


(6.30b) 


for k=l,2,...A. With the <H> given by eq. 6.10, eq. 6.30 become 



(6.31) 


dt 


m 
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(6.32) 


The equations 6.31 and 6.32 are the equations of motion for the 
centroids of the wave packets . 

The right hand side of eq. 6.32 is the net force acting on 

it 

the k^ wave packet and it is found by taking the gradient of eq. 
6.10 with respect to giving 


-3<H> 


r 3 h f 2 i 
= z — — 1 
L i 4 2m 1 ©7 J 


da \ ^1 ’ ‘ ‘ ^ 


[ Z 

L 1 J 


_d<v i j > 


j>i «R. 


( 6 . 33 ) 


As discussed in section 6.2.5, since <v^> depends on © ^ and © ^ 
(eq. 6.23, 6.26) and due to the dependence of © on (R^.Rg ...R^), 
<v. .> becomes non-local in character. Therefore, all the terms 
in the potential energy in eq. 6.10 are required to be 

differentiated with respect to R^ and the second bracket on the 
right hand side of eq. 6.33 does not reduce simply to 
Z- d <v. .>/ 9 R, as in the case of point particles or the 

y i 1J K 

frozen wave packets. 

da eq. 6.33 are obtained by differentiating ©^ in eq. 

6.18 with respect to R k? giving 
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do' i (R 1 ,R 2 , . .R a ) 


d R, 


r i f _L r 

' l b 3 l*i[ a* la. 


i j <?R 


a \ 


+ ^3- 


.2 J 


da 


aR k 


x exp [ - R i / 7a\ ] 


(6.34) 


where, R. . = R. -R . . 

l j l j 


Thus at any given instant t, we first find i 0 '^} from eq. 6.18 
with the instantaneous configuration {R^(t)} and then similarly 
solve the A-coupled eqs . 6 . 34 to find 9a ^ or i = l»2,...A. 
Similarly, we find da for k=l,2,...A. Thus the eqs. 6.31 and 
6.33 together with a ^ and dc^/dR^ as described above and 
<?<v. ,>/dR found from eq. 6.23 and eq. 6.26 completely determine 
the time evolution of a many-particle gaussian wave-packet system 
in this semi-classical microscopic approach . 
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6.3 RESULTS AND DISCUSSIONS 

The equations in section 6.2.6 can be solved in principle to 

find the dynamic evolution of any A-particle system . However , due 

to the iterative solutions of the equations for cr. and do /dR , 

1 1 

and also due to the non-local character of <v, . computations may 

1 3 

become tedious in practice for large values of A and require 

several orders of magnitude more computer time compared to the 

corresponding calculations in the earlier chapters. Therefore, in 

the present work, as an example of the dynamic evolution, we study 

the simplest case of A=2 only. The A=2 case is easier to solve 

and it demonstrates the effects of "Pauli repulsion". Also ? in 

section 6.3.2 we give an example of a many-particle system such as 
1 8 

0 and find its ground state configuration. 

6.3.1 Study of a Two-Body System 

The expectation value <H> (eq. 6.10) for a two wave packet 
system is given by 

<P X > 2 < p 2 >2 3 fi 2 1 3 h 2 1 

<H> = + + p + “ o + < v 1? > (6.35) 

2m 2m 4 2m o-f 4 2m a 

i 7. 

For a two wave packet system it is easy to find an expression for 
the widths o^ and cr ^ of the two wave packets. We note that due to 
the s ymm etry of the system the widths of the two wave packets have 
the same values, i.e. cr - = a. Then^eq. 6.35 reduces to 
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<P 1 > 2 <P 2 > 2 3 h 2 2 


<H> = 


2m 


2m 


+ “ ( —Z ) + <v 1C) > (6.35) 

4 2m o* 1 * 


and eq. 6.18 gives 


o' = bo [ exp (-R 2 2 /2c? 2 )] 1/3 


or ? solving for a we get 


2 

O' = 


or, 


- r: 


12 


2 In b ' 1 


<y - a R 


12 


CG-BS) 


where 


and 


■^12 ~ ^ 2 ^ ~ I <r l‘ > _ 


a = (2 In b 3 ) 1/2 


With the value of b given by eq. 6.19, we find a=0.833. 


We now discuss expression (6.35) for two particular cases 
v 12 = v lg* 


of 
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C a) Case of v . _ ■ 0 
2 c 


For = 0 and a given by eq. 6.36, eq. 6.35 gives, 


<H> = 


<P ± >‘ 


2m 


<P > 

2 3 h 

2m 4 2m 


( a* ) 


2 C 2 
a R 12 


( 6 . 37 ) 


and we find 


-0<H> 


3 ft 4 


A 

R 


12 


. 0 2 ^3 

4 2m a R ^2 


(6.38a) 


and 






3 

4 


A 

-4R 


12 


2m 


a 2 R 3 


12 


(6.38b) 


where R^ 2 is a unit vector along (^- ^ ) ■ 

Substituting the eq. 6.38 in eq. 6.32, we find the equations 
of motion for the centroids R^ and R^ respectively, given by 


m 


d^R. 


dt 


-i-b) 


k 12 


r; 


12 


(6.39a) 


and 


m 


d^. 


dt 


3& 




/N 

R 


12 


12 


(6.39b) 


Thus, we find that even though the two-body interaction 
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between "the two wave-packets is zero ( v i 2 = 0 )> ^ wo wave packets 
are acted upon by a repulsive force given by eq. 6.39. In a 
dynamic evolution the centroids of the two wave packets move apart 
and go to oo as time t-» oo ; while the wave packets spread out with 
their width <y increasing according to eq. 6.36, as R^ 2 increases. 

Therefore, the third term in eq. 6.37 can also be considered 
as a "Pauli repulsive potential" between the two wave-packets. 
This term acts to prevent close encounters of the wave packets, as 
it goes to + oo as R^ 2 + 0. 


C fc>3 Case of v . _ * 0 : 

1 £ 


For v^g * 0, we have 


<H> = 


<P 1 >' 


2m 


<V' 


2m 


3 ^ 


(srl 


4 2m v a 2 R^ 2 


+ <v 12 > 


(6.40) 


N 

with o i = a 2 = a, <v^ 2 > in eq. 6.23 reduces to 


~?2> = 


r o [ 


R, 


C r. 


R 12 R 0 


erf 


( 23^) ] exp ( -jr) (6 41) 


where R Q = (r 2 + 4a 2 R 2 2 ) 1/2 

Differentiating eq. 6.41 with respect to R 12> we get 
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We choose an arbitrary set of parameters (Vq = 3200 MeV, C = 

N 

0.8 fm and rQ = 1.2 fm) for <v i2 > anc ^ the ground state 

configuration for the two wave packet system by the method 
described in section 6.2.5. In the ground state of the/ wave 
packet system with the above potential parameters, we find 
R° 2 =1.23 fm, ^=^2=1-03 fm, <r 2 > 1/2 = 1.89 fm , Cv* 1 .^ = -48.9 

MeV, =29.9 MeV and the total energy E = -19.0 MeV. 

The potential <v^ 2 > ( eq .- 6.41) with the above set of 

parameters is plotted in Figure 6.2 as a function of R^ 2 * ^QM 

and (T qm + <v^ 2 > ) are also plotted in Figure 6.2. From Figure 6.2 

we find that ( t qm +<v i 2 > ^ has a at R i2 = 1,23 which is 

equal to -19.0 MeV. Thus, we see that the determination of the 

ground state by the method described in section 6.2.5 indeed 

produces the minimum of <H> for the two wave packets system. It 

N 

is also to be noticed that the minimum of (Tq^ + <v^ 2 >) is not 

N 

the same value of R 12 as that for <v 12 >. 



Two wave -packet system 
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The dynamic evolution of the two wave packets described 
above, with the ground state as the initial condition, does not 
produce any change in its state. Therefore, the ground state 
obtained above is indeed an equilibrium configuration. 

Starting with the following initial conditions given by <P^> 
= 0, <Pg> = an d Rj 2 “ we study the time development of 

the two wave-packet system. The total energy <H> and the total 

-> -v 

linear momentum (<P^> + <P£>) conserved. Figure 6.3 shows the 
time evolution of the centroids (R^ and the distance R^ 

between the centroids, and the width of the wave packets a . We 
see from Figure 6.3 that R^g oscillates between R^g = 2.0 fm and 
R ^2 “ 0.8 fm corresponding to the respective turning points in 
Figure 6.2. a also oscillates in phase with R^ and with the same 
frequency. 

4 ft 

6.3.2 Study of a Many-Body System C Ground State of CD 

In this section we find the ground state configuration of 

1 6 

0, as an example of the application of the semi-classical method 
described in section 6.2.5. We neglect the spin degree of freedom 
and consider each nucleon as described by a gaussian wave packet 
(eq. 6.4). We choose the NN-potential and the Coulomb potential 
given by expressions (6.23 and 6.26). We also consider the 
strength of the NN interaction between the like nucleons to be 
less than "that between the unlike nucleons by about 20% as 
considered in chapters 4 and 5. 



(fm) 
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Time evolution of the centroid positions R^ and Rg, centroid 
separation R^g and the width and cr ^ of the two 
wave-packets with the potentials shown in the Figure 6.2. 



We choose = 1.2 fro and adjust the parameters Vq and C in 

eq. 6.23 such that the binding energy and the rms radius of 0 in 

its ground state are close to th^corresponding experiments values . 

Parameters of this potential (eq 6.23) and the corresponding 

16 

ground state properties of 0 are given in Table 6.1. 

The distribution n(^) of the width {& i ) of the wave packets 

1 6 

in the ground state of 0 is shown in Figure 6.4 and the radial 
distribution n(r) of the position of the wave packets is given in 
Figure 6.5. The value of c/ ^ varies from 0.78 fm for a wave 
packet in the central core region of the nucleus to about 0.99 fro 
for a wave packet on the surface. The average value of o, is 
about 0.92 fm. 

The calculated value of the average kinetic energy per 

1 6 

nucleon (<T>/A) for the above ground state configuration of 0 
nucleus is about 19 MeV, which is in close agreement with the 
experimentally determined value (16± 3 MeV) from ref. 7. 

The nuclear density distribution is calculated as described 
in chapter 5 section 5.2.4. The total nucleon density (p= p p + p & ) 
is plotted in Figure 6.6. Comparison with the density 
distributions calculated in chapter 5 (Figure 5.4) shows that the 
density distribution claculated in the present case is similar to 

/CJ 

that with the potential <W2> in chapter 5, but it is much smoother 
also. 

Finally, we calculate the compressibility coefficient (K) for 
the ground state of ^^0 . Calculation of K is also similar to 



169 


TABLE 6.1 


Parameters of the effective NN potential and ground-state 
1 B 

properties of 0 



Calculated 

Experimental 

v o £rnev) 

3150.0 


C (fm) 

1.01 


r Q (fm) 

1.20 


Total Nuclear Potential 
energy (MeV) 

-440.2 


Total Coulomb energy (MeV) 

12.9 

7 

<T>/A (MeV) 

18.8 

16± 3 

Binding energy C^eV) 

-127.1 

-127.6 

< r 2 > i/2 

2.73 

2.73 


K (MeV) 


230 



o o 



.7 0.8 0.9 1.0 

<r ( fm) 

FIGURE 6.4 The distribution (n(<?)) of the number of wave-packets 
with widths a in the ground- state of 
potential given in the Table 6.1. 


0 for the 
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r (fm) 

FIGURE 6.5 The distribution (n(r)) of the number of wave-packets 

with their centroids located at radius (r) 

1 6 

ground state of 0. 


in the 
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that, described in chapter 5 section 5.2.4. However, at each 

value of the nuclear radius we calculate the total energy after 

redetermining all the from eq. (6.18) for the compressed or 

the expanded configuration as the case may be. In Figure 6.7 we 

plot the average kinetic energy per nucleon (<T>/A), the average 

potential energy per nucleon V = (Z Z <v. . >/A) and the total 

i j > i 1J 

energy per nucleon (E/A) as a function of the rms radius (R = 
2 1/2 

<r > ) of the compressed or the expanded configuration around 

the ground state. The calculated value of the compressibility 

1 6 

coefficient K for 0 in the present calculation is about 230 MeV. 
This value is slightly higher than the value calculated for 0 in 
chapter 5 section 5.3.2 with the potential <W2> and it is also 
slightly higher compared to the RPA and the constrained 
Hartree-Fock calculations given in Ref. 8. 

In Figure 6.8 we plot Cv!^) 5 Tqm anc * (Tq^ + <v i2 > ^ w ith "the 
parameters given in Table 6.1 for a two wave-packet case 

(deuteron). It is seen from the Figure 6.8 that (Tq^ + <V 12 > ^ 

} 

doesnot have any bound state for a two wave-packet system. 
Therefore, in the present case the deuteron is unbound with the 
potential parameters given in table 6.1. The reason for this 

discrepancy may be in the explicit neglect of the spin degree of 
freedom and the exchange forces which are neglected in the 
present approach. It may also be necessary to consider 
non-spherical gaussian wave packets (eq. 6.1 and 6.2) with crx^, 
ay^ and since the transverse components of the width of two 

wave-packets may not change with respect to the inter-centroid 
separation R^- Thus reducing their contribution to the total 



<r *T <«"> 

FIGURE 6.7 

Plot of the average kinetic energy per nucleon 

(<T>/A), average potential energy per nucleon 

(V/A), and the total energy per nucleon versus rms 

radius (R) of the compressed or expanded 

1 B 

configurations of 0 around R = R 
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Potential <v 12 >, Q.M. kinetic energy and (<v" 2 > + T^) Vs. the centroid 
separation R 12 in the case of a two wave-packet system with the potential 
given in Table 6.1. 
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repulsive kinetic energy (T^). 

6.4 SUMMARY AND CONCLUSIONS 

In the semi-classical microscopic approach presented 
above we assumed that each nucleon in a nucleus can be described 
by a minimum uncertainty gaussian wave packet. We further assumed 
that the wave packets always remain gaussian, but removed some of 
the assumptions made earlier in chapter 5 regarding the width of 
the wave packets. Using the Fermi-gas kinetic energy and the 
local density approximation we determined the width of every 
individual wave packet in the nucleus and its quantum mechanical 
part of the kinetic energy. We found that the width of a wave 
packet depends on its location in the instantaneous configuration. 
Thus, in this approach all the nucleons respect the uncertainity 
relation^ and the effect of the "Pauli repulsion” arising from the 
exclusion principle in a many-body fermionic system is also 
incorporated . 

Assuming the quantum mechanically derived expectation value 
<H> of a many-body system in the Hartree-ansatz, as a new 
classical many-body Hamiltonian, we discussed the dynamic 
evolution of the expectation values of the position {R^} and 
momentum {P^} of the wave packets. 

As a simple example, we studied the dynamic evolution of a 
two wave packet system and found that even in the absence of any 
two body interaction (v 12 = 0), the two wave packets move apart 
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due to the "Pauli-repulsion" - the wave packets spread out as they 
move apart. The zero-point fermi energy or the "Pauli-repulsion*’ 
for the two wave-packet system is found to vary as inverse square 
of the distance between the centroids of the two wave packets. 

In a potential well the centroids of the two-wave packets 
show an oscillatory motion between the two turning points at the 

same total energy. The widths of the wave packets also oscillate 

with the same frequency and in phase with the motion of the 

centroids . 

1 6 

We also determined the ground state configuration of 0 

nucleus by cyclically minimizing the total energy with respect to 

small variations of the individual centroid coordinates. By 

reproducing the ground state binding energy and the rms radius of 
1 6 

0, by adjusting the potential parameters, we find that the 

average kinetic energy per nucleon in the ground state 

1 6 

configuration of 0 is also in close agreement with the 

7 

experimental value. The calculated compressibility coefficient 
has a reasonable value compared to the other theoretical 
calculations and the nuclear density is also smooth. 

It may be worthwhile to mention here that E.J. Heller (Ref. 

and. j.o 

9^ also developed two semiclassical methods for the dynamic 
evolution of the gaussian wave packets. One of the method is 
based on the Ehrenfest theorem which, however, is good only for a 

9 

locally quadratic potential. The other method is based on the 

lo 

Hartree-ansatz and the variational principle. Recently, the 
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second approach was extended by Singer and Smith (Ref. 11) for 
gaussian wave packet simulation of some many-particle molecular 
systems. However, in the approach of Ref. (11) the exchange 
effects and quantum mechanical correlations for a fermionic system 
are neglected. Also, the widths of all the gaussian wave packets 
are assumed to be equal in Ref. 11. Though we also use the 
Hartree-ansatz for the many -body wave functions, the important 
effect of the exclusion principle is incorporated in the present 
approach through the use of a degenerate Fermi-gas kinetic energy 
function . 

• Aichelin and Stocker (Ref. 12) have also used gaussian wave 
packets in a high-energy heavy-ion collision calculation. 
However, in the calculation of Ref. 12 the width of the wave 
packets are assumed to be constant and have the same given value 
for all the wave packets representing individual nucleons. 

The semi-classical microscopic approach described in the 
present chapter for the dynamic evaluation of wave packets can be 
applied in the study of heavy-ion collisions. However, due to the 
enormity of the numerical task, it may take few orders of 
magnitude more computational time compared to the corresponding 
calculations presented in the earlier chapters. As discussed in 
section 6.3.2, it may be necessary to include the spin degree of 
freedom in the interaction potential to reproduce correct nuclear 

€ 

property's for a range of nuclei. Also, it may be necessary to 
consider non-spherical gaussian wave packets, as discussed in 
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section 6.3.2 and also due to the fact that in high energy 
heavy-ion collisions high degree of asymmetry is involved as in 
the case of "side-splash" of nuclear matter (section. 4.3.2)- 
Nevertheless, the semiclassical approach developed here deserves 


further studies. 
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APPENDIX A 


INITIAL COULOMB TRAJECTORIES OF HEAVY-K>NS 


In this appendix we give the formulae for finding the 

positions and velocities of the ions in the centre-of-mass 

reference frame, at a given ion~ion separation (R ) . 

CM 

Let b be the impact parameter for an A 1 +A„ collision, and let 

- 1 * La 

V be the initial velocity of the ion (A.^) with mass moving in 

the ( -x) -direction as shown in the Figure A.l. In Figure A.l, C 

t 

and C are the two focii of the hyperbola that corresponds to the 
Coulomb trajectory of the ion A 1 . The origin of the C.M. frame is 
established at C. 

The components of the position (x^,y^) and velocity 

\ ~ v j*) of ion (A ± ) are given by, 

x 4 = r cos e 
7 ± - r sin 6 

v = v cos e - v sin 6 
\ r y 

Vy = v r sin e + v^cos e 
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vT 



FIGURE A. 1 

Schematic diagram showing the Coulomb trajectory of 
one of the ion in the cent re-of -mass reference 
frame. C and C’ are the two focii of the hyperbola 
that corresponds to the trajectory of the ion. 
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where 


r = R 


CM 


M 


M 


6 = y - y 

o 


Cosy' = 1/e 
o 


Cosy' 


a(£ -1 ) 


£ L 



dr 

dr 

dy 

dt 

dy 

dt 


M 

r 

2 

~ 1 

l 

J 


V b £ Sin y> 


a (* -1) 




dy' 

dt 


M 

-03 


V b 


M A 


reduced mass : M = 


° (M 4 +m 2 ) 


With, 
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r 


2b .2 

J *1 1/2 

Eccentricity : 

' = L 1 

+ ( 

■-) 

] 

Collision parameter: 

a = | 

M 

o 




1 


where, d is the distance of closest-approach for a head-on 
collision (b=0), and it is given by 


d 


ZlZ.ZJe 2 

M V 2 
o 


0 

The components of the position R 2 = (x 2 ,y 2 ) and verity 

v = (v , v ) of the other colliding ion (A,) are given by 

2 x? y2 2 






REFERENCE: 


R.D. Evans, r/ie Atomic Nucleus (McGraw Hill, 1955), Appendix B. 
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APPENDIX B 


EVALUATION OF DOUBLE-FOLDING INTEGRALS FOR THE TWO-BODY POTENTIALS 


<v ij ( l r i~ r jl ]> = i’^i (r i )0 j (r j )v ij ( l r i“ r <j l d3r i d3r j 

.... ( 1 ) 

with (p’s given by 


(r,-<r, >) 




(2,^) 3/4 


exp 


[ 71 + t <P 1 >-(r 1 -<r i »] 


. ( 2 ) 




we get 


1 (r i -<r i >)‘ 

3~ JJ ex P f- 2 1 v iJ ( l r i" r jl ) 

(2no i cy j ) 1 ^i J 3 3 


<v ij> = 


exp 


|r j~ <i j >)i 

a*? 


] d3r i d3r j 


( 3 ) 


Let us, for simplicity, choose the origin of the coordinate system 
at <r i > = 0, and let <rj> = r, then we get 
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<v ij > = 


1 _ r 7 (r,-r) Z 

-3 IS ex P [- ( 2 ~ + 2 ) ] 

v * a *) 1 K 2a f 2o , J J 


( 27X0- ^cr J ) 


v ij (| r 1 - r jl ) d’rjd^j .. (4) 

Now, consider the argument of the exponential in the above 
integral (eq. 4). Multiplying and dividing this term by (o^+ o^) , 
we get 


r 2 (r . -r )‘ 

fi + 

*'1 




[( r i + J 4 ) + [<v r,2+ ;r (r j- r,z )] 


2 (a 


b; IK * If ' 3 


- 2r^ (r j-r ) 


„ 0^0 j 

<r J' r)Z ) + 2r i (r J' r) ~ I 


+ I (r.-r) + 


2(o 


1 r 2 r°^ ^ I 2 1 

2^27 [< r r r i- r > + br r i + b <rj " r) J J 

i J 


(5) 
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Let us make the following transformati 


ons 


0,2 = (cr* + a 2 .)/ 2 


( 6 ) 


R = r j ‘ r i 


(7a) 


R = 


°i 

^ r i + 


(7b) 


d 3 r i d 3 r j = |Jj d 3 R dV 


( 8 ) 


The Jacobian |J| for the above transformations is easily found 


°i°j 


m- (-rV] 
1 J 


(9) 


With the transformations (eq. 7) and substitution (eq.6), the 
identity in eq. 5 becomes 


r\ (r,-r) 2 


2j? 


4 > 

T-- [ 


2 2 
(R-rT R’ ^ 


4 o'' 


4c?‘ 


] 


( 10 ) 


Therefore, substituting eq. 8 and 10 with eq. 9, the double 
folding integral eq. 5 rjluces to a product of two- integrals over 
R’ and R respectively : 


v = Ir ” Ir 


( 11 ) 


where- 
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-R ' 2 

I R , = J exp £ — y ] d 3 R’ 

= (* . 4o 2 ) 3/2 (12) 

and 

r r -(R-r ) 2 t , 

x r = J exp — v( | R| )d*R . 

L 4o- -* 

_ r 2 „ _ r 2 

= exp f — r- 1 J dR R 2 exp f 2~] v{ l R l ) * 

4p * o '•4 o J 


Tl r 

f dft Sin 6 . exp 

o L 


- (R 2 -2Rr Cos©) , 2 n 


4<y 


] J 

J o 




Carrying out the <t> and 6 integrations, we get 


Qncr‘ 


I R 


Let 


r -r -v oo c -R ^ f r 1 

exp [ J Jv(|R|)R exp [— | J sinh l~2 R J dR 


(13) 


= ft 


(14) 


2o f ' 


Substituting eq. 12 and 13 for integrals I R , and I R respectively, 
in eq. 11, and also making use of the substitution given by eq. 

14, we get 



ley 


<v i J > = 


r oo r ^ 

2 J V (|R| )R Sinh(/3R) expf r-1 dR 

v 2o ■'o i 4c ^ J 


. . . (15) 


Now, we evaluate the above integral for the NN potential and 

the Coulomb potential. 

(a) NN Potential 


v i j ( |R| ) = ~ v 0 (1_C/R) ex p (-R 2 /^) 


(16) 


Substituting eq. 16 in eq. 15, we get 


N 

<V ij >: 


Yn a r 


_ r 2 r ^ -R^ ~R^ 

[-^5] [/ R exp [ —2-} exp [^r] 


00 . r ^ 1 f ® 1 

C J Sinh(/5R) exp [ £ J exp [ Tj 


(17) 


( 18 ) 


with this substitution, eq. 16 becomes 
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R 

<v i } > = ;= — «p 

yn a r 


t^f 


co 


J R Sinh (/SR) exp (~a 2 R 2 ) dH 


00 


- c J Sinh (PR) exp (-« 2 R 2 ) 


dR 


(19) 


The two integrals in the above equation can be found in the tables 
of integrals, such as in Ref. 1 or evaluated easily, the results 


are 


00 

J R Sinh (/SR) exp (-a 2 R 2 ) dR = 


4ot‘ 


n r 


.2 6XP l 

4ct 2 J 


( 20 ) 


and 


V* 


Sinh (/SR) exp (-a 2 R 2 ) dR = | — exp [ ~ — £ J erf 


2 x 2 

where erf (x) = — — f e dt, is the error function 

Vn o 


( 21 ). 


'ij* 


integrals 

eq. 20 and 

21 in 

eq. 

19, 

we get 

-V 

o 

2 

r " r i I 

’ M 

71 

exp 

r * 2 1 

= _ exp 

Vno-r 

l 7ff 2 J 1 

■ 4c 2 j 

a 2 

1 4a 2 J 




(22) 



191 


= -V 


■[ 7Z, 


(2aa)' 


r (2 era) 


f 1 

l r ' r 

— 

exp ? + — ~ 

l 2a J 

J 2 a 2 4a ZJ 


(23) 


Substituting for a and ft from eq. 14 and 18, respectively, in eq. 
23, and simplifying the argument of the exponential, we get 



where r = | <r^> " an< ^ a ^- ven ^y ©9- 6. 
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CbD COULOMB POTENTIAL 


r 1.44 

V ij (R) = (MeV) 


(25) 


Substituting in eq. 15, we get, 


<v C > = 

< ij 


Vnor 


exp 


r -r - oo 1.44 . 

( ~zr ) J [— ) B = inh <f3R) 




R 


-R 


exp 


f — 1 dR C^V) 


(26) 


1.44 


-/nor 


^ T) W 

r -r ^ 00 f -K ^ cmeV) 

exp ^ 2 ~ J J Sinh (/3R) ex P ^2j (27) 


4o- ‘ 


From eq. 21, we get, 


<v C > 
<v i j 


1.44 
-/no 'r 


exp 


& { 


2o exp erf ifto) 


) 


C meV} 

(28) 



(29) 


Reference : 

_ j r m Pv^hlk Tables of Integrals , Series, ar«i 

I.S. Gradshteyn and I.M. Ryznix, / 

Products (Academic Press, New York, 1965). 
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